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ABSTRACT 


The  effects  of  an  oscillating  close>coupled  canard  on  the  canard/wing  vortex 
interaction  for  increased  lift  enhancement  were  studied.  Two  test  conditions  were 
studied:  the  first  with  a  model  angle  of  attack  of  22“  and  the  second  of  34“.  The 
canard  was  positioned  at  three  mean  deflection  angles  equal  to  4“,  7°  and  10°  for  the 
model  angle  of  attack  of  22“  and  -4“,  -7“  and  -10“  for  the  model  angle  of  attack  of 
34“.  At  each  of  the  canard  mean  deflection  angles,  the  canard  was  oscillated  with 
amplitudes  of  ±5“  and  ±10“  with  reduced  frequencies  ranging  from  0.046  to  0.232. 
Because  of  the  small  effects  noted  which  were  of  the  order  of  accuracy  of  the 
balance,  only  general  trends  are  discussed.  The  trends  indicate  that  for  this 
particular  model  configuration  and  geometry,  lift  was  decreased  slightly  with 
increasing  canard  frequency  and  amplitude.  No  lift-enhancement  benefits  were 
revealed  during  the  study. 
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I .  INTRODUCTION 


A.  BACKGROUND 

With  the  continued  technological  advances  and  human 
physiological  limitations  becoming  the  edge  of  the  operating 
envelope  in  air  combat,  innovative  technologies  must  be 
developed  in  order  to  maintain  a  decided  advantage  over 
potential  adversaries.  This  advantage  includes  the  ability  to 
engage  an  adversary  and  bring  the  nose  to  bear  upon  the 
adversary  as  quickly  as  possible.  According  to  Dornheim  [Ref. 
1],  [the  development  of]  "All  aspect  weaponry  meant  that  tail 
chase  tactics  no  longer  had  to  dominate  short  range  combat. 
Instead  just  pointing  at  the  adversary  is  sufficient".  As 
HcAtee  [Ref.  2]  states,  "the  decision  to  engage  is  based  on 
time  to  kill  which  has  dropped  considerably  since  World  War 
II".  This  type  of  maneuvering  facilitates  the  use  of  all¬ 
aspect  weapons  where  the  bottom  line  is  being  the  first  one  to 
maneuver  the  aircraft  into  the  weapons  envelope.  Current 
technologies  being  investigated  to  achieve  these  goals  are 
thrust  vectoring  and  post-stall  lift  enhancement  utilizing  a 
close-coupled  canard.  Thrust  vectoring  is  being  studied  on 
the  X-31  fighter  aircraft,  and  the  SAAB  Viggen  is  generally 
credited  as  the  first  operational  aircraft  to  successfully 
Incorporate  a  close-coupled  canard  with  a  delta  wing. 
According  to  Herbst  [Ref.  3],  "The  combination  of  a  delta  wing 
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with  a  properly  designed  and  integrated  close-coupled  canard 
furthermore  improves  maximum  lift  with  less  penalties  of 
canard-wing  interaction.  At  the  same  time,  a  canard  can  be 
used  as  a  control  device  for  optimum  lift-to-drag  throughout 
the  flight  envelope".  Thus,  the  motivation  for  the  study  and 
development  of  post-stall  lift  enhancement  and  control  is 
established. 

Post-stall  lift  enhancement  can  be  seen  with  many 
mechanisms  including  a  close-coupled  canard.  However, 
"dynamic  stall"  as  a  post-stall  lift  enhancer,  which  is  of 
interest  for  this  research,  has  its  roots  in  helicopter  rotor 
dynamics . 

B.  DYNAMIC  STALL 

As  Carr  [Ref.  4]  states,  dynamic  stall  was  first 
encountered  "...when  helicopter  design  engineers  were  unable 
to  predict  the  performance  of  high  speed  helicopters  using 
conventional  aerodynamics".  This  phenomenon  was  also  found  in 
the  flow  over  insect  wings,  compressor  blades  and  the  like. 
When  investigated  it  was  discovered  that  the  rapid  pitching 
experienced  by  the  retreating  blade  caused  the  shedding,  from 
the  leading  edge  (in  accordance  with  the  Biot-Savart  law),  of 
a  strong  dynamic-stall  vortex  which  energized  the  boundary 
layer  causing  the  flow  to  remain  attached  at  angles  of  attack 
well  beyond  the  static  stall  angle  of  attack  for  the  airfoil. 
This  phenomenon  is  transient  in  nature  and  the  lift 
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enhancement  generated  by  the  dynamic-stall  phenomenon 
discontinues  as  soon  as  the  vortex  passes  downstream  of  the 
airfoil.  Carr  [Ref.  4]  studied  dynamic  stall  of  two- 
dimensional  airfoils  in  detail  and  stated  that  "the  need  for 
basic  research  in  three-dimensional  dynamic  stall  effects, 
compressibility  effects  on  dynamic  stall  and  positive  control 
of  unsteady  separated  flow  as  well  as  other  fundamental  areas 
of  unsteady  aerodynamics"  needs  to  be  conducted.  Carr  [Ref. 
4]  also  reviewed  previous  research  on  dynamic  stall  and  found 
that  free-stream  Mach  number  as  well  as  mean  angle  of  attack, 
pitch-axis  location,  amplitude  and  frequency  of  oscillation 
affect  dynamic  stall. 

C.  CANARD  COMFIGDRED  AIRCRAFT 

Lacey  [Ref.  5]  performed  extensive  studies  on  close- 
coupled-canard  location,  sweep  angle  and  size  with  respect  to 
the  main  wing  to  maximize  constructive  interference  and  to 
ultimately  increase  lift.  Results  presented  in  Reference  5 
indicate  that  for  favorable  interference  to  occur,  the  40% 
exposed  root  chord  position  of  the  canard  must  be  placed 
within  1.5  wing  chords  of  the  wing  quarter  chord  with  the 
exposed  trailing  edge  of  the  canard  slightly  forward  of  the 
exposed  leading  edge  of  the  wing  root.  The  canard  must  be 
placed  above  the  wing  with  an  optimum  vertical  separation  of 
0.1  to  0.25  wing  chords.  Canard  leading  edge  sweep  should 
meet  or  exceed  60°.  The  area  ratio  (canard  area/wing  area  = 
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(S^/S^))  that  provided  the  maximum  lift  coefficient,  = 
0.25,  was  also  the  largest  value  tested,  implying  studies  of 
larger  area  ratios  may  need  to  be  conducted.  Lift  was 
decreased  for  small  angles  of  attack  in  all  cases,  but  lift 
was  increased  beyond  the  static  stall  angle  of  the  wing  alone 
at  higher  angles  of  attack. 

Further  research  conducted  by  Er-El  [Ref.  6]  supported 
the  conclusions  of  Lacey  [Ref.  5]  and  added  insight  into  the 
dynamics  of  the  canard-wing  interaction.  Er-El  [Ref.  6] 
stated,  "The  changes  in  the  leading-edge-vortex  trajectories 
due  to  the  presence  of  the  canard  are  evident  in  the  pressure 
distributions.  Comparison  of  the  spanwise  position  of  the 
suction  peaks  shows  that  in  the  wing/canard  configuration 
these  peaks  are,  in  general,  outboard  to  those  in  the  wing 
alone  configuration".  Er-El  [Ref.  6]  also  commented, 
"Evidently,  the  interaction  between  the  wing  and  canard 
vortices  in  the  aft  section  of  the  wing  is  stronger  when  the 
canard  is  highly  swept  due  to  its  stronger  vortices" . 

D.  PREVIOUS  STUDIES 

Because  dynamic  stall  is  transitory  in  nature,  creative 
engineering  is  necessary  in  order  to  exploit  the  benefits  of 
dynamic  stall  on  an  as-needed  and  continuing  basis.  Because 
dynamic  stall  is  created  by  the  shedding  of  a  strong  leading- 
edge  vortex  following  the  rapid  pitch  of  an  airfoil,  the 
airfoil  must  continually  be  pitched  or  oscillated  to  provide 


4 


the 


the  desired  lift  enhancement, 
oscillation  requirement, 


Due  to  this  pitch  or 
close-coupled-canard 


configuration  is  useful  because  the  canard  can  be  pitched  in 
an  oscillatory  manner  continuously  to  shed  the  leading-edge 


vortex  necessary  to  possibly  energize  and  reattach  the  flow 


over  the  main  wing  providing  post-stall  lift  enhancement. 


Interest  has  been  generated  in  determining  whether  a  highly 


swept  canard  located  optimally  as  determined  by  Lacey  [Ref.  5] 
and  oscillated  continuously  would  help  reattach  the  flow  over 
the  main  wing,  thus  utilizing  the  dynamic  stall  phenomenon  to 
enhance  post-stall  maneuverability. 

Ashworth  et  al.  [Ref.  7]  obtained  qualitative  and 
quantitative  data  when  the  canard  of  an  X- 2 9 -type  ha If -body 
aircraft  was  sinusoidally  oscillated.  Smoke  visualization 
showed  the  canard  leading-edge  vortex  passing  well  above  the 
main  wing  and  the  canard  tip  vortex  being  split  by  the  main 
wing.  Hot-wire  measurements  above  and  below  the  wing 
indicated  that  the  flow  velocity  above  the  wing  between  20% 
and  40%  chord  was  slightly  less  than  the  free-stream 
velocity,  indicating  the  existence  of  a  separation  bubble. 
Below  the  wing  the  velocity  was  essentially  the  free-stream 
value.  No  lift  and  drag  measurements  comparing  the  results  of 
a  non-oscillating  canard  and  the  results  of  an  oscillating 
canard  were  taken. 


5 


Huyer  et  al.  [Ref.  8]  conducted  similar  studies  where 
flow  visualization  and  surface-mounted  pressure  transducers 
were  used  to  characterize  the  flow  fields  over  the  trailing 
airfoil.  This  information  was  used  to  compute  the  pressure 
distributions  as  well  as  the  normal  and  tangential  force 
coefficients  with  the  trailing  airfoil  at  angles  of  attack  of 
10°  and  20°.  Huyer  et  al.  [Ref.  8]  analyzed  the  pressure 
distribution  data  and  discussed  the  unsteady  nature  of  the 
pressure  distribution.  Additionally,  Huyer  et  al.  [Ref.  8] 
discussed  the  canard  dynamic-stall  vortex  effects  on 
canard/wing  interaction  and  stated  "The  production  of  the 
[dynamic  stall]  vortex  did  however,  appear  to  energize  the 
boundary  layer  and  redirect  the  potential  flow  over  the 
trailing  airfoil.  This  phenomenon  resulted  in  flow 
reattachment  over  the  trailing  airfoil  surface  at  angles  of 
attack  far  exceeding  static  stall  angles."  Huyer  et  al. 
[Ref.  8]  however,  does  not  compare  the  oscillating  canard  data 
to  fixed  canard  data  nor  is  there  any  lift  or  drag  information 
presented . 

Hebbar  et  al.  [Ref.  9]  utilized  flow  visualization  to 
compare  wing  vortex  burst  location  of  an  X-31-like  fighter 
aircraft  with  a  fixed  and  oscillating  close-coupled  canard. 
Reference  9  concludes  that  small  amplitude  and  low  frequency 
oscillations  cause  early  vortex  bursting  whereas  the  small 
amplitude  high  frequency  oscillations  appear  to  delay  vortex 
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bursting.  It  is  worth  noting  that  the  reduced  frequencies 
studied  in  Reference  9  (lc=1.7  to  k==10.4),  are  well  above  the 
reduced  frequencies  associated  with  dynamic  stall  (k=0.1  to 
k=0.3)  studied  in  this  research. 

To  this  point  previous  research  has  focused  on  the 
qualitative  aspects  as  well  as  the  transitory  effects  of 
canard  oscillation  on  canard/wing  vortex  interaction.  Mabey 
et  al.  [Ref.  10]  addresses  an  important  point  when  he  asks 
"can  dynamic  movement  of  the  canard  (about  a  constant  mean 
setting)  alter  the  mean  flow  on  the  canard  and  hence  alter  the 
mean  flow  on  the  wing?”  This  question  is  the  focus  of  this 
research  and  quantitative  measurements  of  lift  and  drag  were 
taken. 

B.  STATEMENT  OF  PURPOSE 

The  purpose  of  this  research  was  to  build  on  previous 
research  done  by  Kersh  [Ref.  11]  and  Schmidt  [Ref.  12]  in 
considering  the  canard/wing  interaction.  Kersh  [Ref.  11] 
studied  the  baseline  lift  and  drag  characteristics  of  a  low- 
aspect-ratio  wing/body  model  from  low  angles  of  attack  to 
angles  of  attack  beyond  36°.  Additionally,  Kersh  [Ref.  11] 
studied  the  comparative  lift  enhancement  using  the  same 
wing/body  and  a  close-coupled  canard  for  various  wing/body 
angles  of  attack  and  various  canard  deflection  angles. 

The  results  of  Kersh 's  [Ref.  11]  studies  were  used  to 
determine  the  test  conditions  studied  by  Schmidt  [Ref.  12] 
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based  on  the  maximum  lift  enhancement  obtained  with  the  close- 
coupled  canard  configuration.  Schmidt  [Ref.  12]  obtained  lift 
and  drag  measurements  comparing  the  difference  between  tests 
in  which  a  fixed  canard  and  an  oscillating  canard  were  used  to 
determine  whether  significant  increases  in  lift  could  be 
obtained  due  to  the  dynamic-stall  vortex  effects  on  the 
canard/wing  vortex  interaction.  Schmidt  [Ref.  12]  had 
significant  difficulties  with  his  experimental  setup  and  the 
canard-oscillating  mechanism  and  was  not  able  to  examine  the 
full  range  of  test  conditions  desired.  This  research  followed 
the  recommendations  of  Schmidt  [Ref.  12]  for  model 
improvement.  Additionally,  increased  frequency  and  amplitude 
ranges  were  examined  to  complete  the  study  of  the  effects  of 
an  oscillating  close-coupled  canard  on  post-stall  lift 
enhancement . 


8 


II.  EXPERIMENT 


A.  WIND  TUNNEL 

The  Naval  Postgraduate  School  (NPS)  low-speed  wind  tunnel 
located  in  the  basement  of  Halligan  Hall,  as  shown  in  Figure 
1,  was  utilized  for  this  experiment.  As  discussed  in  the  NPS 
Laboratory  Manual  for  Low  Speed  Wind  Tunnel  Testing  [Ref. 
13  ] ,  the  wind  tunnel  is  64  feet  in  length  and  varies  in  width 
between  21.5  and  25.5  feet.  It  is  powered  by  a  100-horsepower 
electric  motor  connected  to  a  four-speed  International  truck 
transmission.  The  transmission  connects  to  a  three-blade 
variable-pitch  fan  capable  of  providing  test  section 
velocities  approaching  200  miles  per  hour.  Directly 
downstream  of  the  fan  are  eight  flow  straighteners  to  reduce 
txirbulence  and  swirl.  The  flow  passes  through  turning  vanes 
at  each  of  the  four  90°  turns  that  complete  the  wind  tunnel 
circuit.  In  the  settling  chamber,  two  fine-wire  turbulence 
screens  six  inches  apart  are  installed  to  help  provide  smooth 
test  section  flows.  The  test  section  has  a  cross  sectional 
area  of  approximately  10  square  feet  (45  inches  wide  and  32 
inches  high)  and  is  approximately  four  feet  long.  The 
contraction  ratio  between  the  settling  chamber  and  the  test- 
section  is  10:1.  Test-section  lighting  and  a  reflection  plane 
located  in  the  test  section  reduce  the  usable  cross-sectional 
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area  to  approximately  9.88  square  feet.  The  turbulence 
intensity  in  the  test  section  is  about  0.2%  and  a  1/500  tunnel 
diameter  breather  slot  is  installed  just  downstream  of  the 
test  section  for  pressure  recovery  due  to  tunnel  leakage. 
Temperature  is  measured  using  a  dial  thermometer  with  the 
probe  located  in  the  settling  chamber.  Test  section  dynamic 
pressure  (g)  is  determined  using  a  water-filled  micromanometer 
whose  readout  is  in  centimeters  of  water.  The  settling 
chamber  and  the  test-section  each  have  a  set  of  four  wall- 
mounted  static  pressure  taps  with  each  set  of  four  connected 
to  individual  common  manifolds  and  then  to  the  micromanometer. 
The  resolution  of  the  micromanometer  was  approximately  equal 
to  ±0.1  cm.  H2O  and  the  micromanometer  could  be  zeroed  to 
compensate  for  evaporation.  Additionally,  a  digital  pressure 
transducer  was  used  for  backup.  The  AP  reading  was  converted 
to  dynamic  pressure  using  a  calibration  equation  from  Kersh 
[Ref.  11]  as  shown  in  Equation  (1).  Test  section  velocity  was 
calculated  using  Bernoulli's  equation  and  the  definition  of 
dynamic  pressure.  Equation  (2) . 

g=2. 047  *( -0.026749+1. 1149AP)  (1) 


1/2+p 


(2) 
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with  the  variables  defined  as  follows: 
p  =  Density  (slugs/ft^) 

AP  =  Micronanometer  reading  (cm.  H2O) 
q  =  Dynamic  pressure  (q  =  27.33  Ibf/ft^) 

V  =  Velocity  (ft/s) 

1.1149  -  Tunnel  calibration  factor 
2.047  =  Conversion  factor  for  cm.  HjO  to  Ibf/ft^ 
-0.02675  =  Tunnel  calibration  intercept 

The  wind  tunnel  calibration  factor  and  tunnel  calibration 
intercept  correct  micromanometer  readings  of  AP  in  the  range 
tested  (3-12  cm.  HjO)  to  the  actual  test-section  dynamic 
pressure.  During  calibration  the  actual  test-section  dynamic 
pressure  was  measured  by  a  pitot-static  probe  place  in  the 
test-section.  The  results  were  plotted  against  the 
micromanometer  readings,  and  with  the  assumption  of  a  linear 
relationship  between  the  pitot-static  probe  and  micromanometer 
readings,  a  slope  (tunnel  calibration  factor)  and  a  Y- 
intercept  (tunnel  calibration  intercept)  were  calculated  using 
a  linear  regression  analysis  as  discussed  in  Yuan  [Ref.  14] 
and  Kersh  [Ref.  11]. 

B.  STRAIN  GAGE  BALANCE  AND  TURNTABLE 

As  discussed  in  Schmidt  [Ref.  12],  the  external  strain- 
gage  balance  was  originally  built  by  NPS  personnel  in  1974. 
It  was  designed  to  measure  normal  and  axial  forces  as  well  as 
normal  and  axial  moments.  The  balance  was  capable  of 
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Figure  2:  Strain  Gage  Balance 


Figur*  3:  Balance  Base  and  KarXings 
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measuring  forces  up  to  150  Ibf  and  each  of  the  four  strain- 
gage  bridge-circuits  had  four  active  legs  with  automatic 
temperature  compensation.  The  forces  and  moments  were 
measured  by  a  pair  of  orthogonally-mounted  strain-gages 
cemented  on  the  flexure  links  of  the  balance  column  and  were 
separated  by  a  vertical  distance  of  26.5  inches  (see  Figure 
2)  .  The  balance  and  turntable  was  capable  of  rotation  from 
-18°  to  +200°  relative  to  the  tunnel  centerline  and  rotation 
was  controlled  remotely  by  an  electric  motor  (see  Figure  3) . 
The  base  of  the  balance  was  marked  in  degrees  so  that  the 
operator  could  determine  the  angle  of  attack  of  the  balance  in 
relation  to  the  tunnel  centerline  (see  Figure  3)  .  When  the 
tunnel  was  in  operation  the  forces  and  moments  generated  by 
the  model  translated  into  voltage  readings  on  each  of  the  four 
strain-gage  bridge-circuits.  These  voltages  were  sent  to  the 
data  acquisition  system  (discussed  in  section  D.  of  this 
chapter)  and  the  normal  and  axial  forces  were  resolved  into 
lift  and  drag  and  stored  on  floppy  disc.  The  model  (discussed 
in  Appendix  A.)  was  mounted  to  the  reflection  plane  with  1/8- 
inch  spacers  to  prevent  the  model  from  rubbing  on  the 
reflection  plane  when  the  angle  of  attack  was  changed  (see 
Figure  4).  The  turntable  had  a  1/8-inch  gap  between  itself 
and  the  reflection  plane  to  allow  the  balance  to  deflect  under 
loading  (see  Figure  4) .  The  balance  was  calibrated  before  use 
to  ensure  accurate  measurements  as  well  as  to  verify  that  the 
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Figure  4:  Turntable  and  Model 


strain-gage  bridge-circuits  responded  linearly  to  loading. 
The  calibration  process  is  discussed  in  detail  in  Appendix  B. 

C.  CANARD/WING  MODEL 

The  canard/wing  model  (see  Figure  5)  was  a  half -body  made 
from  mahogany  and  aluminum  and  mounted  to  an  aluminum  plate 
which  mounted  directly  to  the  balance  turntable  (see  Figure 
4) .  The  model  was  a  generic  fighter-type  aircraft  with  a 
close-coupled  canard  and  low-aspect-ratio  wing  (see  Figure  5)  . 
The  model  consisted  of  three  parts  including  the  ogive  nose, 
which  was  permanently  mounted  to  the  aluminum  base.  The 
second  section  was  the  canard-body  section  which  housed  the 


canard  oscillation  mechanism  and  motor.  On  the  top  of  the 
second  section  where  the  canard  was  mounted  was  a  scale  so 
that  the  canard  mean  deflection  angle  with  respect  to  the 
fuselage  centerline  as  well  as  the  canard  amplitude  could  be 
determined.  The  third  section  or  aft-section  contained  the 
wing  which  was  bolted  to  the  aft  body  from  the  inside  and 
enclosed  the  pick-up  for  the  frequency  counter  used  to 
determine  the  canard-oscillation  frequency.  The  model  design 
is  described  in  detail  in  Appendix  A. 

D.  DATA  ACQUISITION  HARDWARE  AMD  SOFTWARE 

Each  strain-gage  bridge  had  an  individual  signal- 
conditioner  that  provided  the  bridge  excitation  voltage  (+10.0 
volts)  and  a  span  zero  (0.0  volts)  so  that  the  bridge  could  be 
balanced  before  use.  The  output  from  each  bridge  circuit  was 
sent  to  a  Pacific*  8256  low-noise  amplifier  with  the  gain  set 
at  1000.  The  outputs  from  the  amplifiers  were  sent  to  and 
processed  by  a  National  Instriiments*  MC-MIO-16L-9  12-bit 
multi-function  board.  The  MC-MIO-16L-9  has  an  analog-to- 
digital  converter  with  a  9.0  fjtsec.  conversion  time  capable  of 
100  Kbytes/sec.  The  board  had  a  digitation  span  of  4096  bits 
resulting  in  a  4.88-mvolt  resolution  as  discussed  in  the  MC- 
MIO-16L-9  user  manual  (Ref  15].  The  outputs  of  the  MC-MIO- 
16L-9  board  were  then  sent  to  an  IBM  PS-2*  computer.  A  data- 
acquisition  program  written  by  Schmidt  [Ref.  12]  was  modified 
and  utilized  for  automatic  data  acquisition  and  reduction. 
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The  program  was  written  using  Microsoft  QuickBasic  4  •  5* 
implementing  the  MC-MIO-16L-9  board  commands  and  compiled 
using  National  Instruments  LabWindows*.  The  software  was 
written  so  that  each  channel  would  be  scanned  1000  times  over 
a  2.25-second  interval  (444.44  Hz).  The  complete  listing  of 
the  code  and  subroutines  utilized  can  be  found  in  Appendix  C. 
E.  TEST  CONDITIONS 

The  purpose  of  this  experiment  was  to  determine  whether 
an  oscillating  close-coupled  canard  would  provide  significant 
increases  in  lift  when  compared  to  the  results  of  a  static 
close-coupled  canard  with  the  same  configuration.  The  model 
configuration  was  based  upon  research  done  by  Lacey  [Ref.  5] 
in  the  early  1970 's  and  the  test  conditions  to  be  studied  were 
based  upon  the  results  of  Kersh  [Ref.  11]  and  Schmidt  [Ref. 
12].  The  results  of  Kersh  [Ref.  11]  indicate  that  for  a  model 
angle  of  attack  of  22°,  lift  was  maximized  for  a  canard 
deflection  angle  of  +7°  and  for  a  model  angle  of  attack  of  34°, 
lift  was  maximized  for  a  canard  deflection  angle  of  -7°. 
These  maximum  values  were  bracketed  for  this  experiment  by 
adding  or  subtracting  3°  to  the  canard  deflection  angle. 

This  research  follows  directly  work  done  by  Schmidt  [Ref. 
12]  and  Incorporates  model  improvements  discussed  in  Appendix 
A.  The  pivot-point  for  the  canard  was  moved  to  25%  of  the 
mean  aerodynamic  chord  of  the  canard  to  reduce  changes  in 
pitching  moment  with  canard  deflection  angle  thus  avoiding 
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problems  encountered  by  Schmidt  [Ref.  12].  For  all  the  test 
conditions,  a  AP  equal  to  12  cm.  H2O  was  utilized  which 
provided  a  dynamic  pressure  (q)  of  27.33  Ibf/ft^. 

The  first  model  angle  of  attack  (AOA)  studied  was  +22°, 
with  the  canard  at  mean  deflection  angles  equal  to  +4°,  +7°  and 
+10“  with  respect  to  the  fuselage  centerline.  For  each  of 
these  three  test  conditions,  readings  were  taken  with  the 
canard  fixed  and  with  it  oscillating  at  two  amplitudes  ,  ±5“ 
and  ±10“,  at  frequencies  ranging  from  5  hertz  to  25  hertz  in 
increments  of  5  hertz.  The  results  are  presented  in  reduced 
frequencies  and  lift  coefficients  addressed  in  Chapter  IV 
(Discussion  of  Results) . 

The  second  model  AOA  studied  was  +34°,  with  the  canard  at 
mean  deflection  angles  equal  to  -4°,  -7°  and  -10°  with  respect 
to  the  fuselage  centerline.  For  each  of  these  three  test 
conditions,  readings  were  taken  with  the  canard  fixed  and  with 
it  oscillating  at  two  amplitudes,  ±5°  and  ±10°,  with 
frequencies  ranging  from  5  hertz  to  25  hertz  in  increments  of 
5  hertz.  The  results  are  presented  in  reduced  frequencies  and 
lift  coefficients  addressed  in  Chapter  IV  (Discussion  of 
Results) . 
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Ill .  EXPERIMENTAL  PROCEDURES 
A.  PRE-RUN  CALIBRATION 

Before  any  experiments  were  performed,  the  balance  was 
calibrated  and  a  calibration  matrix  was  determined  for  use  in 
the  data  acquisition  program  to  resolve  the  strain-gage  bridge 
outputs  into  lift  and  drag.  The  calibration  procedure  is 
discussed  in  detail  in  Appendix  B  and  reviewed  here  briefly. 


Figurs  6:  cable  and  Pulley  Setup 


The  turntable  was  aligned  with  the  wind  tunnel  visually 
and  the  calibration  rig  was  attached  directly  to  the 
turntable.  Two  cable  heights  above  the  reflection  plane  of 
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Figure  7t  Calibration  Rig 
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7.75  and  10.75  inches  were  utilized  to  resolve  axial  and 


normal  moments.  The  cable  was  strung  over  a  pulley  and 
leveled  with  a  sight-level  and  then  attached  to  a  weight-cage 
(see  Figure  6  and  Figure  7).  Prior  to  beginning  the 
calibration  process  the  inputs  to  the  Pacific*  8255  amplifiers 
were  shorted.  With  the  gain  set  to  1  the  output  screw  was  set 
to  0.0  ±  100.0  ^volts.  The  gain  was  then  set  to  1000  and  the 
input  screw  was  set  to  0.0  ±  500.0  /xvolts.  The  shorts  were 
then  removed  and  the  bridge-circuit  excitation  voltage  was  set 
to  10.0  ±  0.05  volts  and  the  span  voltage  set  to  0.0  ±  0.05 
volts  at  which  point  the  calibration  process  was  started. 
Weights  were  suspended  in  various  increments  up  to  45  Ibf  from 
the  cable  to  apply  a  known  load  to  the  balance  and  turntable. 
This  was  repeated  for  both  cable  heights;  then  the  turntable 
was  rotated  90**  and  the  process  repeated.  The  results  were 
utilized  to  obtain  the  calibration  matrix  used  during  the  test 
runs.  The  calibration  matrix  enabled  the  data  acquisition 
program  to  resolve  the  strain-gage  bridge-circuit  outputs 
directly  into  lift  and  drag  measurements.  Because  the  span 
voltages  were  very  difficult  to  zero  and  the  amplifiers  were 
prone  to  noise  and  drift  at  a  gain  of  1000  (see  Error 
Analysis,  Chapter  V.),  a  tare  reading  was  taken  before  each 
measurement  and  subtracted  from  the  measurement.  The  software 
was  written  so  that  each  channel  would  be  scanned  at  444.44 
hertz  for  2.25  seconds  with  the  readings  averaged,  resulting 
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in  a  single  value  displayed  and  recorded  to  floppy  disc.  Each 
point  was  measured  four  times  to  ensure  consistency  and 
repeatability.  The  response  from  the  strain-gage  bridge- 
circuits  proved  to  be  linear  as  expected,  and  the  calibration 
matrix  obtained  compares  closely  with  that  obtained  by  Schmidt 
[Ref.  12].  The  calibration  data  can  be  found  in  Appendix  B. 
B.  DATA  COLLECTION 

The  test  conditions  as  discussed  in  Chapter  II,  Section 
E,  were  utilized  and  are  repeated  for  clarity: 

1.  Measured  AP  -  12  cm.  H2O 

2.  Test  section  velocity  *=  150  ft/s  (Approximately) 

3.  Model  AOA  =22° 

Canard  Mean  AOA  =  +4°,  +7°,  +10° 

Amplitude  =  ±5°,  ±10° 

Frequency  =0-25  hertz 

4.  Model  AOA  =34° 

Canard  Mean  AOA  =  -4°,  -7°,  -10° 

Amplitude  =  ±5°,  ±10° 

Frequency  =0-25  hertz 

Since  the  canard/wing  model  had  been  significantly 
redesigned  (see  Appendix  A.),  no  baseline  validation  of 
Schmidt's  [Ref.  12]  data  was  attempted. 

A  BK  Precision*  voltage  supply  with  variable  voltage  and 
current  settings  was  utilized  to  provide  power  to  the  canard- 
oscillation  motor.  Canard-oscillation  frequency  was  measured 
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using  a  Monsanto*  Model  lOlA  universal  counter  which  counted 
the  teeth  of  a  gear  mounted  on  the  motor  drive  shaft.  One 
revolution  of  the  motor  drive  shaft  was  equal  to  one  full- 
cycle  of  canard  oscillation  and  by  dividing  the  number 
displayed  on  the  counter  by  the  number  of  teeth  on  the 
sprocket,  a  frequency  in  hertz  was  obtained. 

Because  of  difficulties  encountered  by  Schmidt  [Ref.  12], 
a  set  screw  was  installed  in  the  model  to  maintain  the  mean 
deflection  angle  for  the  static  test  runs.  After  data  were 
taken  at  the  static  test  condition,  the  tunnel  had  to  be  shut 
down  and  the  set  screw  removed  before  the  canard  was 
oscillated. 

A  pre-start  checklist  was  utilized  for  each  test 
condition  and  is  as  follows: 

1.  Operational  amplifier  input/output  zeroed 

2.  Excitation  voltage  set  and  span  voltage  zeroed 
Note:  The  first  two  procedures  were  completed  only 
once  before  each  block  of  tunnel  runs  (i.e.  daily). 

3.  Tunnel  temperature  recorded. 

4.  Atmospheric  pressure  recorded. 

5.  Micromanometer  zeroed. 

6.  Model  configuration  set  and  checked  (with  set  screw) . 

7.  Tare  reading  taken. 

8.  Tunnel  brought  up  to  12  cm.  H2O. 

9.  Fixed-canard  data  point  measured  four  times. 
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10.  Tunnel  shut  down  and  set  screw  removed. 


11.  Tunnel  brought  up  to  12  cm.  HgO. 

12.  Oscillatlng-canard  data  points  measured  four  times 
each. 

13.  Zero  reading  taken  after  tunnel  shut  down. 

14.  Model  reconfigured  and  steps  3-14  repeated. 

The  model  and  canard-oscillating  mechanism  performed 
extremely  well  (see  Acknowledgments)  and  no  difficulties  were 
encountered  in  setting  the  test  conditions  desired  or 
collecting  data  at  the  prescribed  test  conditions. 
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IV.  DISCUSSION  07  RESULTS 


The  following  chapter  covers  all  the  test  cases  as 
discussed  previously,  and  graphical  results  are  presented  in 
reduced  frequencies  which  were  calculated  using  Equation  1. 

(1) 

2*V 

where : 

e  s  Canard  oscillation  frequency  (radians/sec) 

=  5.38  in.  (0.4483  ft.) 

V  =  Test  section  velocity  (ft/s) 

The  lift  and  drag  coefficient  calculations  were  performed 
using  Equations  2  and  3  below. 

r-lift  (2) 


^  _  drag  (3) 

q*S 

where: 

lift  =  lift  (Ibf) 
drag  «  drag  (Ibf) 

S  »  Wing  reference  area  to  the  fuselage  centerline 
(135.12  in^  or  0.9383  ft^) 

The  reduced  experimental  data  can  be  seen  in  Appendix  D. 
The  reduced  data  show  the  mean  values  at  each  test  condition 
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as  well  as  the  standard  deviation  for  each  test  condition. 


The  mean  values  are  the  data  used  for  the  graphical 
presentations . 

A.  MODEL  AOA  =  22°,  AMPLITUDE  s  15° 


OSCILLATING  CANARD 

MODEL  kOk--22,  AMPLITUDE  =  5 


WING  AREA  TO  FUSELAGE  CENTERLINE 
USED  FOR  Srel  IN  Cl  CALCULATION 

Pigure  8:  Model  AOA  =  22  deg. 

Figure  8  shows  the  final  results  for  the  canard  mean 
deflection  angle  equal  to  +4°,  +7°  and  +10°  with  an  amplitude 
of  ±  5°.  An  initial  examination  of  the  reduced  data  showed 
the  curve  for  a  canard  mean  deflection  angle  of  +7°  indicating 
an  increase  in  the  lift  coefficient  at  5  hertz  (k-0.046)  when 
compared  to  the  static  case.  Thereafter,  as  the  canard 
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REPEATABILITY  EXPERIMENT 

MODEL  AOA  =  22;  CANARD  AOA  =  7 


— “  data  set  l  — data  set  2  data  set  3 

data  set  4  data  set  5  — original  data  set 

WING  AREA  TO  FUSELAGE  CENTERLINE 
USED  FOR  Sref  IN  Cl  CALCULATION 

rigur*  9t  Repeatability  Experiaent 

oscillation  frequency  increased  beyond  5  hertz  (ksO.046)  the 
lift  coefficient  decreased.  Because  of  the  uniqueness  of  this 
increase,  the  5  hertz  (k=0.046)  data  point  was  of  interest; 
repeatability  experiments  were  attempted  without  success. 
The  repeatability  experiment  was  conducted  in  the  same  manner 
as  was  the  original  experiment.  For  each  test  run  the  set 
screw  was  installed  and  data  taken,  then  the  tunnel  was  shut 
down  and  the  set  screw  removed  for  the  dynamic  run.  Figure  9 
shows  the  results  of  the  repeatability  experiments  where  all 
but  the  original  data  reflect  a  decrease  in  the  lift 


29 


coefficient  at  5  hertz  (k=0.046).  The  hypotheses  for  the 
discrepancy  in  results  of  the  repeatability  experiments  lie  in 
the  nature  of  the  experimental  setup.  As  discussed 
previously,  the  static  test  had  to  be  conducted  with  a  set 
screw  in  the  model  to  hold  the  canard  mean  deflection  angle. 
After  the  static  run  the  tunnel  was  shut  down  and  the  set 
screw  removed.  The  tunnel  was  restarted  and  the  dynamic 
readings  were  taken.  Because  the  tunnel  had  to  be  shut  down 
between  the  static  and  dynamic  runs,  there  was  an  opportunity 
for  a  bias  error  to  be  introduced.  The  micromanometer 
utilized  has  a  readability  of  approximately  ±  0.1  cm.  H2O 
which  equated  to  a  AC^  equal  to  0.01  (discussed  in  Error 
Analysis,  Chapter  V.).  While  Ac^  equal  to  0.01  is  not  a  large 
change,  it  is  on  the  same  order  of  magnitude  as  the  increase 
in  the  lift  coefficient  observed.  Each  data  point  had  a 
statistical  mean  and  standard  deviation  introducing  other 
errors  which  are  also  discussed  in  detail  in  the  Error 
Analysis  chapter. 

Because  of  the  results  obtained  from  the  repeatability 
experiments,  the  decision  was  made  to  adjust  the  second  data 
point  at  5  hertz  (k=0.046),  as  well  as  all  subsequent  data 
points  from  10  to  25  hertz  (k=0.092  -  0.23)  downward  by  a  AC^ 
equal  to  0.039.  The  data  in  Appendix  D  reflect  the  original 
values  obtained  and  not  the  corrected  values;  the  values 
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displayed  in  Figure  8  were  corrected  for  presentation  and 
discussion  purposes. 

with  the  correction  applied  for  the  canard  mean 
deflection  angle  equal  to  7°,  the  test  results  indicate  that 
the  lift  coefficient  decreases  steadily  with  an  increase  in 
canard  frequency  and  the  trend  is  consistent  for  all  three 
test  cases. 

B.  MODEL  AOA  =  22",  AMPLITUDE  =  ±10” 

Figure  10  shows  the  test  results  for  canard  mean 


OSCILLATING  CANARD 

MODEL  AOA=22,  AMPLITUDE=10 


WING  AREA  TO  FUSELAGE  CENTERLINE 
USED  FOR  Srei  IN  Cl  CALCULATION 


Figure  10 t  Model  AOA  =  22  deg. 
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deflection  angles  equal  to  +4°,  +7“  and  +10“  with  an  amplitude 
of  ±10°.  The  same  general  trends  as  seen  in  Figure  8  are 
present  in  that  the  lift  coefficient  decreased  with  canard 
frequency  regardless  of  canard  mean  deflection  angle.  It  is 
interesting  to  note  that  the  decrease  in  lift  coefficient  is 
more  dramatic  between  the  static  and  the  first  dynamic  test 
condition,  5  hertz  (k=0.046),  when  compared  to  the  results  of 
the  ±5°  amplitude  case  (Figure  8) .  The  decrease  in  lift 
coefficient  is  consistent  for  all  three  test  conditions. 

C.  MODEL  AOA  =  34°,  AMPLITUDE  s  ±5° 

Figure  11  shows  the  test  results  for  canard  mean 
deflection  angles  equal  to  -4°,  -7°  and  -10°  with  an  amplitude 
of  ±5°.  The  lift  coefficient  decreased  significantly  between 
the  static  and  first  dynamic  case,  5  hertz  ()c=0.046),  with  the 
canard  mean  deflection  angle  equal  to  -7°  and  -10°,  then  shows 
a  relatively  flat  response  to  frequency  for  the  remainder  of 
the  frequency  range  tested.  For  a  canard  mean  deflection 
angle  equal  to  -4°,  Figure  11  shows  an  initial  increase  up  to 
about  10  hertz  (k=0.092)  then  a  decrease  with  increased 
frequency  for  the  remainder  of  the  frequency  tested. 

D.  MODEL  AOA  =  34°,  AMPLITUDE  s  ±10° 

Figure  12  shows  the  test  results  for  canard  mean 
deflection  angles  equal  to  -4°,  -7°  and  -10°  with  an  amplitude 
of  ±10°.  Similarly  to  the  ±5°  amplitude  case,  the  canard  mean 
deflection  angles  equal  to  -7°  and  -10°  show  a  significant  drop 
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OSCILLATING  CANARD 

MODEL  AOA=34,  AMPLITUDE  =  5 


WING  AREA  TO  FUSELAGE  CENTERLINE 
USED  FOR  Srel  IN  Cl  CALCULATION 

rlgur*  lit  Mod«l  aoa  s  34  dag. 

between  the  static  and  first  dynamic  test  case  of  5  hertz 
(k«0.046).  The  results  then  show  a  relatively  flat  response 
throughout  the  remainder  of  the  frequency  range  tested.  For 
a  mean  canard  deflection  angle  equal  to  -4°,  the  response  is 
relatively  flat  throughout  the  frequency  range  tested. 

B.  SUMMARY  OF  RESULTS 

For  the  model  AOA  equal  to  22°,  oscillating  the  canard 
seemed  to  have  a  negative  effect  on  the  lift  coefficient. 
However,  as  is  discussed  in  the  Error  Analysis  chapter,  the 
margin  for  error  in  the  measurements  is  as  large  as  the 
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OSCILLATING  CANARD 

MODEL  AOA^S'I,  AMPLITUDEMO 


WING  AREA  TO  FUSELAGE  CENTERLINE 
USED  FOR  Sref  IN  Cl  CALCULATION 

Figur«  12s  Model  AOA  s  34  dog. 

relative  difference  in  the  measurements.  This  prohibits 
precise  interpretation  of  the  results,  but  the  general  trends 
seem  to  indicate  that  the  lift  coefficient  is  adversely 
affected  when  the  canard  is  oscillated,  and  the  canard  dynamic 
stall  vortex  negatively  affects  the  canard/wing  vortex 
interaction  and  these  negative  effects  are  increased  with 
Increasing  frequency.  The  increased  canard  amplitude  of  ±10° 


seemed  to  decrease  the  lift  coefficient  further  because  of  the 
larger  drop  (between  the  static  case  and  first  dynamic  case) 
than  for  the  ±5°  amplitude  case.  For  the  model  AOA  equal  to 


34°,  the  lift  coefficient  is  adversely  affected  the  most  for 
canard  mean  deflection  angles  equal  to  -7°  and  -10°.  The  trend 
is  relatively  flat  after  the  first  dynamic  case  of  5  hertz 
(k-0.046),  indicating  relative  frequency  independence.  For 
the  canard  mean  deflection  angle  equal  to  -4°,  a  slight 
increase  in  lift  coefficient  for  the  amplitude  equal  to  ±5°  is 
seen  up  to  10  hertz  (k=0.092),  then  a  decrease  for  the 
remainder  of  the  frequency  range  tested  is  noted.  For  the 
amplitude  of  ±10°  there  is  a  flat  response  for  the  entire 
frequency  range  tested.  Care  must  be  taken  in  interpreting 
the  increase  in  lift  coefficient  for  the  canard  mean 
deflection  angle  equal  to  -4°  with  the  amplitude  equal  to  ±5° 
because  the  increase  is  well  within  the  statistical  error 
limits  as  discussed  in  the  Error  Analysis  chapter. 


V.  ERROR  ANALYSIS 


Because  of  the  small  magnitude  of  change  between  the 
measurements  taken,  particular  attention  had  to  be  paid  to 
error  analysis.  There  appears  to  be  three  sources  of  error 
that  will  be  discussed  separately.  The  first  source  of  error 
is  the  change  in  g  (Ibf/ft^)  due  to  the  experimental  setup 
requiring  the  tunnel  be  shut  down  between  the  static  and 
dynamic  test  runs.  The  second  source  of  error  is  in  the  data 
acquisition  system  and  the  third  source  is  the  statistical 
distribution  of  the  data.  Finally,  no  consideration  of  a 
fixed  error  or  bias  was  considered  because  information  on  the 
relative  differences  between  data  points  was  sought,  not 
absolute  measurements. 

A.  ONCERTAINTY  IN  DYNAMIC  PRESSURE  (q) 

A  test  was  conducted  to  determine  the  tunnel  sensitivity 
to  lift  coefficient  variations  with  respect  to  changes  in  AP 
(dACi/dAP) .  Figure  13  shows  a  linear  relationship  between 
error  in  lift  coefficient  with  respect  to  AP  (cm.  HjO)  .  There 
is  a  positive  slope  because  the  lift  coefficient  calculation 
is  based  on  the  assumption  that  AP,  which  equals  12  cm.  H2O, 
is  constant  which  is  consistent  with  the  data  reduction 
procedures.  This  linear  relationship  has  a  positive  slope 
with  a  value  of  dAC^/dAP  equal  to  0.1185.  When  the 
readability  of  the  micromanometer  is  considered  (AP  =  ±0.1  cm 
H2O) ,  it  is  evident  that  if  for  the  static  test  case  the 
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TUNNEL  SENSITIVITY 

TO  DELTA  IN  PRESSURE 


DELTA  P  (cm  H20| 


Cl  =  .1185  *  cm.  H20 
actual  data 

MODEL  AOA  =  22;  CANARD  AOA  =  2 
12  cm.  H20  USED  AS  REFERENCE 

Figure  13: Tunnel  Sensitivity 

micromanoineter  were  at  the  lower  end  of  this  interval  and  for 
the  first  dynamic  case  the  micromanometer  reading  were  at  the 
high  end  of  this  interval  an  error  in  the  lift  coefficient 
(AC^)  as  great  as  0.02  is  possible.  This  case  is  an  extreme 
one  and  is  unlikely  but  must  be  considered.  Because  the 
tunnel  was  not  shut  down  between  the  dynamic  runs  (5  -  25 
hertz,  k>»0.046  -  0.23)  this  error  will  not  be  present  in  the 
dynamic  readings.  Specifically,  the  relative  difference 
between  the  dynamic  readings  would  be  accurate  within  the 
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statistical  and  data  acquisition  limitations  and  not  subject 
to  differences  in  q  (Ibf/ft^). 

B.  STATISTICAL  SIGNIFICANCE  OF  DATA 

Because  of  the  small  sample  size  of  the  data  taken  at 
each  test  condition,  special  statistical  procedures  were 
utilized  to  statistically  legitimize  the  data.  For  a  normal 
distribution  with  a  large  sample  size  the  central  limit 
theorem  as  discussed  in  McClave  [Ref.  16]  states,  "If  a  random 
sample  of  "n"  observations  is  selected  from  a  population  (any 
population) ,  then,  when  "n”  is  sufficiently  large,  the 
sampling  distribution  of  will  be  approximately  a  normal 
distribution.  The  larger  the  sample  size,  "n",  the  better 
will  be  the  normal  approximation  to  the  sampling  distribution 
of  This  theorem  lays  the  foundation  necessary  for  the 

calculation  of  the  mean  and  standard  deviation  from  which 
confidence  intervals  for  the  results  can  be  established.  But 
when  presented  with  small  sample  sizes,  work  done  by  W.S. 
Cosset  as  discussed  in  McClave  [Ref.  16],  states  that  if  a 
random  sample  is  selected  from  a  population  with  a  normal 
distribution,  the  sampling  distribution  will  be  approximately 
normal.  However,  the  variance  of  this  small  sample  will  have 
a  greater  uncertainty  associated  with  it,  so  the  t-statistic 
was  developed.  It  accounts  for  larger  uncertainty  in  both  the 
mean  and  the  variance  due  to  the  small  sample  size,  so  when 
developing  confidence  intervals  they  are  much  larger  than  for 
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the  normal  z-statistic  as  discussed  in  McClave  [Ref.  16].  The 
t-statistic  uses  an  "n-l”  sample  size  and  for  a  95%  confidence 
interval,  the  standard  deviation  is  multiplied  by  3.182  which 
compares  to  1.96  for  the  z-statistic.  The  upper  limit  and  the 
lower  limit  for  the  confidence  interval  as  well  as  the 
variance  can  be  found  in  the  reduced  experimental  data  in 
Appendix  D.  Confidence  intervals  as  large  as  ±0.03  can  be 
found  which  are  on  the  same  order  of  magnitude  as  the  changes 
in  the  lift  coefficient  observed. 

C.  DATA  ACQUI8ITI0H  UNCERTAINTY 

Because  of  problems  encountered  by  Kersh  [Ref.  11]  and 
Schmidt  [Ref.  12],  the  data  acquisition  system  was  examined 
for  its  resolution  capability.  Additionally,  it  was  known 
that  the  data  acquisition  outputs  tended  to  drift  cyclically 
with  time.  Figure  14  shows  how  the  value  of  the  lift 
coefficient  drifted  over  a  period  of  1/2  hour.  The  tunnel  was 
started  and  brought  up  to  a  AP  equal  to  12  cm.  H2O  and  data 
were  taken  every  five  minutes  for  30  minutes.  The  AP  data  was 
plotted  versus  time  and  a  linear  regression  was  performed. 
The  variance  between  the  predicted  value  obtained  from  the 
regression  equation  and  the  actual  value  at  a  specific  time 
were  calculated  to  specify  a  confidence  interval.  The 
assumption  that  the  drift  in  the  system  is  linear  for  short 
periods  of  time  can  be  seen  to  be  reasonably  accurate  with  a 
90%  confidence  interval.  The  drift  of  the  acquisition  system 
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TUNNEL  DRIFT 

AS  A  FUNCTION  OF  TIME 


— ^  lower  limit  — linear  regression 

upper  limit  a  actual  data 


90  %  CONFIDENCE  INTERVAL 

Figur*  14:  Tuim«l  Drift 

was  determined  to  be  Ac^  equal  to  -0.03  ±0.01  per  hour  for 
this  case.  As  a  note,  because  the  drift  is  oscillatory  in 
nature,  the  drift-curve  slope  could  just  as  easily  have  been 
positive  and  quite  probably  would  have  been  if  the  drift  were 
observed  at  some  other  time.  The  primary  purpose  of  this 
discussion  is  to  establish  that  the  data  acquisition  system 
does  in  fact  drift  over  the  duration  of  a  typical  test  run. 

Because  the  drift  is  not  linear  and  is  more  oscillatory 
over  large  periods  of  time,  only  a  30-minute  period  was 
examined  because  each  test  run  performed  between  tare  readings 
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took  approximately  30  minutes.  The  source  of  drift  was 
thought  to  be  from  amplifier  output- voltage  drift.  This  drift 
was  not  applied  to  any  of  the  results  discussed  in  Chapter  IV. 

Because  of  the  unpredictable  nature  of  the  drift  as  well 
as  the  fact  that  the  first  type  of  error  can  be  applied  only 
between  the  first  two  data  points,  a  standard  uncertainty 
analysis  as  discussed  in  Holman  and  Gajda  [Ref.  17]  cannot  be 
conducted.  The  importance  of  this  discussion  is  evident  when 
quantitative  conclusions  about  the  data  are  to  be  reached. 
Because  the  uncertainty  is  on  the  same  small  order  of 
magnitude  as  the  changes  in  lift  coefficients,  definite 
quantitative  comparisons  in  the  results  cannot  be  made. 
However,  general  trends  in  the  results  can  be  noted  to  the 
effect  that  oscillating  the  canard  provided  no  significant 
increase  in  lift  as  expected. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


An  oscillating  close-coupled  canard  was  studied  to 
determine  the  effect  on  the  canard/ wing  vortex  interaction  for 
increased  lift  enhancement.  Two  test  conditions  were  studied: 
the  first  with  a  model  angle  of  attack  equal  to  22°  and  the 
second  equal  to  34°.  The  canard  was  positioned  at  three  mean 
deflection  angles  equal  to  4°,  7°  and  10°  for  the  model  angle 
of  attack  equal  to  22°  and  -4°,  -7°  and  -10°  for  the  model 
angle  of  attack  equal  to  34°.  At  each  of  the  canard  mean 
deflection  angles,  the  canard  was  oscillated  with  an  amplitude 
equal  to  ±5°  and  ±10°  with  reduced  frequencies,  k,  ranging  from 
0.046  to  0.232.  The  following  conclusions  were  drawn  from  the 
test  results. 

-  Because  of  the  small  effects  noted  which  were  of  the 
order  of  accuracy  of  the  balance,  only  general  trends  can  be 
noted . 

The  trends  indicate  that  for  this  particular  model 
configuration  and  geometry,  lift  was  decreased  slightly  with 
increasing  canard  frequency  and  amplitude. 

The  following  recommendations  are  made: 

-  Further  studies  of  different  model  geometries  with  an 
oscillating  canard  need  to  be  conducted.  These  studies  should 
include  varying  the  pivot-point  of  the  canard,  varying  the 
dynamic  pressure  (q) ,  and  varying  wing  and  canard  geometries 
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as  well  as  varying  the  frequency  and  amplitude  of  canard 
oscillation. 

-  The  canard-oscillation  mechanism  needs  to  modified  so 
that  the  mean  angle  of  attack  can  be  maintained  without  a  set 
screw.  This  will  eliminate  the  necessity  of  shutting  the 
tunnel  down  between  the  static  and  dynamic  test  runs,  thus 
eliminating  one  source  of  error. 

-  The  strain-gage  balance  and  data  acquisition  system 
need  to  be  replaced  with  a  system  capable  of  more  accurate 
measurements.  The  balance  itself  is  old  and  has  been  abused 
in  the  past,  and  the  outputs  of  the  bridge  circuits  will  be  in 
question  as  the  balance  is  continually  used. 

-  Flow  visualization  could  be  used  to  further  understand 
the  canard/ wing  vortex  interaction. 
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APPENDIX  A:  MODEL  DESIGN 


The  model  characteristics  and  design  parameters  were 
based  on  work  previously  done  by  Schmidt  [Ref.  12],  Kersh 
[Ref;  11]  and  Lacey  [Ref.  5]  and  incorporated  recommendations 
by  Schmidt  [Ref.  12]  for  the  canard-oscillation  mechanism. 


A.  GENERAL  MODEL  CONSIDERATIONS 


The  model  can  be  seen  in  Figures  A1  through  A3  and 
utilizes  the  same  fuselage  as  Kersh  [Ref.  11]  and  Schmidt 
[Ref.  12].  An  NACA  64A008  airfoil  section  was  used  for  both 
the  main  wing  and  the  canard.  No  attempt  to  trip  the  boundary 
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Figure  A2:  Model  Configuration 

existing  aircraft  airfoil  shapes.  The  model  was  mounted  to 
the  balance  turntable  17.18  inches  back  from  the  model  tip  of 
the  ogive  nose,  and  the  model  was  36  inches  long.  The 
fuselage  was  4.5  inches  wide  and  3  inches  high  with  a  semi¬ 
span  of  14.6  inches  as  measured  from  the  reflection  plane  to 
the  wing  tip.  The  canard  trailing  edge  did  not  overlap  the 
wing  as  discussed  in  Lacey  [Ref.  5]  and  the  canard  mean 
deflection  angle  was  monitored  by  lines  drawn  on  the  model  at 
the  trailing  edge  of  the  canard  in  one-degree  increments  with 
the  zero-degree  deflection  angle  corresponding  to  the  fuselage 
centerline. 
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Figure  A3:  Model  Configuration 

the  zero-degree  deflection  angle  corresponding  to  the  fuselage 
centerline. 

B.  CANARD  MODIFICATIONS 

The  canard  size  and  shape  were  determined  by  Kersh  [Ref. 
11]  based  upon  work  done  by  Lacey  [Ref.  5]  and  Behrbohm  [Ref. 
18]  and  included  a  leading-edge  sweep  equal  to  60“,  taper 
ratio  (A)  equal  to  0.1  and  an  aspect  ratio  (AR)  equal  to  2  as 
calculated  by  Equations  A1  and  A2. 


AR=2 


C^d+A.) 


^  (A2) 

where: 

b  =  Wingspan 

=  Exposed  root  chord 
=  Tip  chord 

The  canard-oscillation  mechanism  was  redesigned  and 
reconfigured  to  handle  the  dynamic  loadings  created  by  canard 
oscillations  as  discussed  by  Schmidt  [Ref.  12],  and  the 


Figure  A4:  canard  Oscillation-Mechanism 
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mechanism  can  be  seen  in  Figure  A4.  The  redesign  allowed  for 
changes  to  the  mean  deflection  angle,  amplitude  and  frequency 
of  the  canard  oscillation.  The  mean  deflection  angle  was  set 
by  adjusting  the  turnbuckle  attached  to  the  canard  shaft-arm 
and  the  eccentric-hub.  The  amplitude  was  adjusted  by  changing 
the  hole  in  which  the  turnbuckle  was  connected  in  the 
eccentric-hub,  and  fine  adjustments  were  made  by  sliding  the 
turnbuckle  up  or  down  in  the  canard  shaft-arm.  Frequency  was 
controlled  by  the  motor  and  a  variable  power  supply.  The 
pivot-point  for  the  canard  was  moved  from  7%  MAC  (40%  exposed 
root  chord)  to  25%  MAC  to  reduce  changes  in  pitching-moment 
with  changes  in  deflection  angle  and  thus  reduce  the  loading 
on  the  oscillation  mechanism.  Also,  past  work  has  involved 
dynamic-stall  configurations  pivoted  about  25%  MAC;  this 
modification  would  put  the  current  study  in  line  with  previous 
efforts.’  This  change  moved  the  canard  forward  on  the  model 
and  changed  the  basic  configuration  determined  to  be  optimal 
by  Lacey  [Ref.  5]  with  x/c  less  than  or  equal  to  1.5,  where  x 
is  the  distance  from  40%  of  the  canard  exposed  root  chord  to 
25%  MAC  of  the  wing  and  c  is  the  MAC  of  the  main  wing  as 
determined  by  Equation  A3. 

UJ) 
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To  compensate  for  the  new  configuration,  a  new  wing  was 
designed. 

C.  WING  MODIFICATIONS 

The  wing  maintained  the  same  basic  dimensions  as  the  wing 
used  by  Schmidt  [Ref.  12],  with  a  leading-edge  sweep  equal  to 
50°,  taper  ratio  (1)  equal  to  0.15  and  an  aspect  ratio  (AR) 
equal  to  3,  but  the  overall  size  was  changed  from  68.89  square 
inches  to  90.19  square  inches  (exposed  area)  as  seen  in  Figure 
A5.  This  larger  main  wing  resulted  in  an  x/c  equal  to  1.44 
which  was  within  parameters  as  specified  by  Lacey  [Ref.  5]. 
Additionally,  previous  research  done  by  Lacey  [Ref.  5],  Kersh 
[Ref.  11]  and  Schmidt  [Ref.  12]  used  main-wing  areas  to  the 
fuselage  centerline.  Because  the  fuselage  is  very  large  in 
order  to  house  the  canard-oscillation  mechanism,  distorted 
area  ratios  resulted  because  the  bulk  of  the  projected  wing- 
area  was  inside  the  fuselage.  Therefore,  the  exposed  wing- 
area  was  used  when  sizing  the  airfoil  and  resulted  in  a  ratio 
of  the  canard  exposed-area  to  the  wing  exposed-area  (S^^/S^) 
of  0.21.  The  vertical  separation  between  the  canard  and  wing, 
z/c,  was  equal  to  0.21  where  z  was  the  vertical  distance 
between  the  main  wing  and  the  canard.  This  parameter  is 
within  the  optimum  range  for  z/c  of  O.l  to  0.25  as  determined 
by  Lacey  [Ref.  5]. 
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APPENDIX  B:  BALANCE  CALIBRATION 


Procedures  as  outlined  in  Schmidt  [Ref.  12]  were  utilized 
to  obtain  the  calibration  matrix  to  be  utilized  in  the  data 
acquisition  program  (see  Data  Acquisition,  Appendix  C) . 
Additionally,  it  was  sought  to  duplicate  the  calibration 
operation  and  to  obtain  a  similar  calibration  matrix,  thus 
validating  the  results  obtained  by  Schmidt  [Ref.  12]. 

An  externally-mounted  cylindrical-type  balance  with 
strain-gages  located  at  the  base  (pair  A)  and  the  top  (pair  B) 
of  the  balance  has  been  discussed  previously.  Strain-gage 
pair  A  and  B  were  separated  by  a  vertical  distance  of  26.5 
inches  and  temperature  compensation  was  accomplished  utilizing 
four  active  legs  for  each  strain-gage  bridge.  A  schematic  of 
the  balance  and  strain  gage  location  can  be  seen  in  Figure  Bl. 
The  top  of  the  balance  or  txirntable  was  flush  with  the 
reflection-plane  in  the  test-section  of  the  low-speed  wind 
tunnel  and  a  calibration  rig  was  mounted  to  the  turntable  for 
the  calibration  process  as  shown  in  Figure  Bl.  With  the 
Pacific*  8255  operational  amplifiers  set  to  a  gain  of  1000, 
and  the  MC-MIO-16L-9  board  set  to  a  gain  of  1,  the  system  was 
capable  of  analog-to-digital  conversion  with  a  4.88-mvolt 
resolution.  With  the  calibration  rig  installed,  the  angle-of- 
attack  readout  at  the  base  of  the  balance  was  visually  set  to 
ensure  that  the  angle  of  attack  read  on  the  base  corresponded 
to  the  actual  model  angle  of  attack.  The  inputs  to  the 
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Figure  Bl:  Strain-Oage  Balance  Scbematic 
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operational  amplifiers  were  shorted  and  the  gain  set  to  1;  the 
amplifier  output  was  then  adjusted  to  0.0  ±  100.0  /xvolts.  The 
gain  was  then  set  to  1000  and  the  amplifier  input  was  adjusted 
to  0.0  ±  500.0  /xvolts.  The  shorts  were  removed  and  the  signal 
conditioners  were  set  to  a  span  voltage  of  10  ±  0.05  volts  and 
the  span  zeroes  were  set  to  0.0  ±  0.05  volts. 

Because  the  balance  was  only  capable  of  rotation  from 
-18°  to  +200°,  and  the  installation  of  the  model  required  the 
turntable  to  rotate  from  o“  to  -90°,  the  turntable  was  rotated 
90°  so  that  a  90°  reading  on  the  turntable  base  (see  Figure  B2) 
corresponded  to  a  model  angle  of  attack  equal  to  0°.  Because 
of  this  rotation,  the  normal  and  axial  force  measurements  were 
rotated  90°  as  seen  in  Figure  B3.  The  data  acquisition  and 
force  reduction  equations  accounted  for  this  offset  where  the 
90°  position  corresponded  to  the  0°  position  in  the  new 
coordinate  system.  The  nomenclature  for  each  strain-gage  is 
as  follows: 

Eaa  -  Voltage  at  the  lower  axial-force  bridge. 

Eba  -  Voltage  at  the  upper  axial-force  bridge. 

Ean  -  Voltage  at  the  lower  normal-force  bridge. 

Ebn  -  Voltage  at  the  upper  normal-force  bridge. 

(a  -  b)  -  Height  above  turntable  of  first  cable 

attachment  (hs=10 . 7 5  in . ) . 

(a*  -  b)  -  Height  above  turntable  of  second  cable 

attachment  (h*7.75  in.). 
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Figure  B2:  Balance  Base  and  Markings 


Figure  B3:  Coordinate  System  cbange 


57 


The  calibration  procedure  consisted  of  rotating  the 
turntable  to  either  0°  or  90°,  as  measured  in  the  new 
coordinate  system,  and  suspending  increasing  increments  of 


Figxire  B4:  Calibration  Rig  and  Pulley 


weights  on  the  calibration  rig  (see  Figure  B4  and  B5)  at  two 
heights  above  the  turntable  and  using  the  data  acquisition 
system  to  record  the  voltage  readings  obtained.  The  cable 
height  was  measured  vertically  from  the  reflection  plane  to 
the  cable  and  a  sight-level  was  used  to  ensure  the  cable  was 
horizontal. 
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Figure  B5:  Calibration  Rig  and  Weight  Cage 
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Equation  B1  is  the  basic  equation  used  to  determine  the 
calibration-matrix ; 


[tC] 
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Where  [K]  and  [dA/dLOAD]  are  4X4  matrices.  The  1X2  sub¬ 
matrices  on  the  right  hand  side  are  for  moment  resolution  and 
the  entire  equation  can  be  expanded  as  follows: 


[dAE„ 

dAE^ 

d^E,, 

dAEi^ 

dA 

dA 

dN 

dN 

■^11 

^13 

^4 

dAEi^ 

dAEy^ 

dAEt^ 

K22 

^3  ^4 

dA 

dA 

dN 

dN 

^2 

^33 

^34 

dAB„ 

dAEia 

dA£„ 

dAE'^ 

[^41 

^43 

^44, 

dA 

dA 

dN 

dN 

dAEj. 

dAE^a 

dAEto 

dAEjtu 

dA 

dA 

dN 

dN 

1 

1 

0 

0 

ia-b) 

ia'-b) 

0 

0 

0 

0 

1 

1 

0 

0 

ia-b) 

ia'-b). 

60 


The  right  side  of  the  equation  is  known  with  (a-b)  equal  to 
7.75  inches  and  (a'-b)  equal  to  10.75  inches.  The  slopes  of 
the  voltage  readings  with  respect  to  load,  cable  height  and 
turntable  angle  of  attack  (0°,  90°)  were  determined  using  the 
calibration  data  and  performing  a  linear-regression  analysis 
to  fill  in  the  slope  matrix  ( [dAE/dLOAD] ) .  Both  sides  of  the 
Equation  B2  were  pre-mult iplied  by  the  inverse  of  the  slope 
matrix  resulting  in  the  determination  of  the  calibration 
matrix  [K] .  The  response  of  the  strain-gages  were  linear 
within  the  load  ranges  tested  and  graphical  results  of  the 
calibration  process  can  be  seen  in  Figures  B6  through  B9. 

The  normal  and  axial  forces  as  well  as  the  normal  and 
axial  moments  were  resolved  using  Equation  B3: 


where : 

=  Axial  force  (Ibf) 

=  Axial  moment  (ft  Ibf) 

Nf  =  Normal  force  (Ibf) 
n„  =  Normal  moment  (ft  Ibf) 

[K]=  Calibration-matrix 

E,,  (etc.)  =  Strain-gage  voltage  output 
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Lift  and  drag  were  calculated  using  Equations  B4  and  B5: 


L=Af*cos  (a)  +Nf*sinia)  (b4) 


D=Af*sin(a)  -Nf*cos(a)  (B5) 


where: 

o  =  model  angle  of  attack 

and  all  other  variables  were  previously  defined  in  Equation 
B3.  The  results  were  recorded  to  floppy  disc. 

The  non-dimensional  axial  and  normal  coefficients  were 
calculated  using  equations  similar  to  Equations  2  and  3  in 
Chapter  IV. 
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Normal  Voltage  Readings 

Average  Calibration  Data 


II  :  7.75  In. 


Figure  B7:  Normal  Calibration  Data 
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Figure  B8:  Axial  Calibration  Data 


Normal  Voltage  Readings 

Average  Calibration  Data 


MIMvolts  (mV)  (Ttousaids) 


(0  15  20  25  30  35 

Weight  (Ibf) 

■Eaa  Eba  -•-Eaa  -"“Ebn 


a  >  10. 7S  («. 


Figure  B9:  Horaal  Calibration  Data 


REDUCED  CALIBRATION  DATA 


AXIAL  VOLTAGE  READINGS  WITH  CABLE  HEIGHT  =  10.75  in. 


1 - 

weight 

Eaa 

Eba 

Ean 

Ebn 

0 

0.002442 

0.200196 

0.50293 

0.313721 

3.204 

488.3191 

196.5137 

5.997314 

4.641113 

8.21 

1298.835 

535.1221 

12.44751 

10.61401 

13.206 

2074.525 

876.8555 

21.7395 

24.01978 

23.209 

3673.109 

1598.728 

32.00317 

54.35669 

33.213 

5339.808 

2454 

42.54395 

168.8855 

36.217 

5703.39 

2650.382 

46.03272 

199.7754 

33.213 

5328.344 

2468.192 

42.15332 

184.7668 

23.209 

3834.091 

1675.983 

36.04126 

57.15942 

13.206 

2110.815 

894.7705 

21.92627 

22.13501 

8.21 

1291.871 

533.2349 

12.18994 

11.13892 

3.204 

500.1721 

203.2727 

6.73584 

4.321289 

0 

0.643311 

-0.16602 

-0.09644 

0.067139 
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REDUCED  CALIBRATION  DATA 


REDUCED  CALIBRATION  DATA 


NORMAL  VOLTAGE  READINGS  WITH  CABLE  HEIGHT  =  10.75  in. 


weight 

Eaa 

Eba 

Ean 

Ebn 

0 

-0.0354 

-0.04639 

0.252685 

0.45166 

1  3.204 

-6.35742 

-15.2698 

508.6426 

342.2449  | 

8.21 

-18.2312 

-33.2935 

1339.437 

900.6689  1 

13.206 

-29.1589 

-28.2129 

2116.339 

1441.707  1 

23.209 

-52.865 

31.71997 

3804.698 

2682.972  1 

33.213 

-73.8831 

16.37329 

5319.387 

3721.821  1 

36.217 

-83.1555 

7.347412 

5872.511 

4095.486  1 

33.213 

-76.9397 

8.337403 

5331.328 

3717.717  1 

23.209 

-54.6875 

-44.6008 

3750.99 

2579.695  | 

13.206 

-30.8728 

-53.4937 

2105.709 

1430.215 

03 

« 

to 

H 

-19.1748 

-35.1257 

1297.524 

877.9834 

3.204 

-7.03491 

-13.6133 

511.5002 

345.2881 

0 

-0.03784 

0.062256 

0.05127 

-0.14648 
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REDUCED  CALIBRATION  DATA 


NORMAL  VOLTAGE  READINGS  WITH  CABLE  HEIGHT  =7.75  in. 


weight 

Eaa 

Eba 

Ean 

Ebn 

0 

-0.08179 

0.146485 

-0.24536 

0.039062 

3.204 

-10.1282 

-12.3987 

475.0464 

276.8579  1 

8.21 

-27.5647 

-37.3169 

1223.66 

711.1499 

18.212 

-62.0129 

-69.071 

2687.106 

1570.194  1 

28.215 

-99.5422 

-62.6099 

4225.94 

2509.735  1 

38.219 

-134.57 

-14.4495 

5684.46 

3433.84  1 

43.229 

-153.827 

-11.1719 

6529.763 

3948.425  I 

47.236 

-165.776 

44.88281 

7072.443 

4319.435  1 

43.229 

-152.38 

-23.4387 

6509.513 

3926.559  I 

38.219 

-136.27 

-44.5374 

5851.049 

3511.797  1 

28.215 

-102.057 

-69.3469 

4316.512 

2562.994 

18.212 

-63.7305 

-73.8782 

2767.905 

1625.876 

8.21 

-28.5706 

-37.1631 

1266.998 

740.8362 

3.204 

-9.98535 

-13.5437 

464.7095 

272.749 

0 

0.024414 

0.2771 

0.050049 

0.91919 
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APPENDIX  C:  DATA  ACQUISITION  PROGRAM 


This  program  was  written  and  compiled  using  LabWindows  and 
QuickBasic  4.5.  (used  "be  /o  multi"  to  compile)  Its 
purpose  is  to  read  and  convert  voltages  from  four  channels 
connected  to  the  strain  gauges  on  the  Academic  wind  tunnel. 
The  voltages  are  converted  to  normal  and  axial  forces  and 
moments  with  respect  to  the  balance.  It  was  written  and 
modified  by  LT  Tom  D.  Stuart  and  LT  Dean  c.  Schmidt,  20  June 
92.  It  was  further  modified  by  LT  Douglas  G.  Me  Bane  on  12 
January  1993.  It  follows  the  nomenclature  used  by  Kersh  where 
the  table  is  rotated  90  degrees  and  the  axial  force  is  normal 
to  the  flow  direction  and  the  normal  force  is  parallel  to  the 
flow  direction  (opposite  direction) .  However,  the  lift  and 
drag  components  are  converted  so  that  the  lift  is 
perpendicular  to  the  flow  direction  and  the  drag  is  parallel 
to  the  flow  direction  (same  direction)  concurrent  with 
standard  notation. 


'  Variables  explained 

'  eaa  =  Strain  gauge  voltage  at  point  A  in  Axial  direction. 

•  eba  =  Strain  gauge  voltage  at  point  B  in  Axial  direction. 

'  ean  =  Strain  gauge  voltage  at  point  A  in  Normal  direction. 

'  ebn  =  Strain  gauge  voltage  at  point  B  in  Normal  direction. 

'  AX  -  Axial  force 

'  Max  ==  Axial  moment 

'  NORM  =  Normal  force 

'  Mnorm  =  Normal  moment 

'  alpha  =  Angle  of  Attack  of  the  model 

'  LIFT  *  Lift  force 

'  DRAG  =  Drag  force 


REM 

REM 

REM 

REM 

REM 

REM 

REM 


$INCLUDE: 
$ INCLUDE: 
$INCLUDE; 
$INCLUDE: 
$INCLUDE; 
$INCLUDE: 
$INCLUDE: 


'  C : \LW\ INCLUDE\LWSYSTEM . INC • 
'C: \LW\INCLUDE\GPIB.INC' 

'  C : \LW\INCLUDE\FORMATIO . INC ' 
'  C : \LW\ INCLUDE\GRAPHICS . INC ' 
•C: \LW\INCLUDE\ANALYSIS.INC' 
•  C : \LW\ INCLUDE\DATAACQ . INC ' 

'  C : \LW\INCLUDE\RS232 . INC ' 


DIM  K#(4,4) 

DIM  ean. array# (1000) , eaa. array# (1000) , ebn. array# (1000) , 
eba. array# ( 1000) 


69 


E 


D 


COMMON  SHAR 

ean .  array#  ( ) ,  eaa .  array#  ( ) ,  ebn .  array#  ( ) ,  eba .  array#  ( ) 

DECLARE  SUB  volt  (ean#, eaa#, ebn#, eba#, alpha#) 

DECLARE  SUB  aero  (AX#, NORM#, LIFT#, DRAG#, alpha#) 
DECLARE  SUB  forces 
(K#  ( )  ,  eaa# ,  eba# ,  ean# ,  ebn# ,  AX# , Max# , NORM# , Mnorm# ,  alpha#) 

SCREEN  9,  0 
COLOR  15,  1 

eaaO#  «  o 
eaaO#  =  0 
eanO#  =  0 
ebnO#  =  0 


*  Ifk  It -kit -kit  It  It -kit  hit  hit  hie  hit  1e  It -kit  It  hit -kit  hit  hie  It -kit -kick  it  It -kit  *1(11  ****** 

•  CALIBRATION  MATRIX  INPUT  (See  thesis  for  explanation) 


DATA  0.009292, 
DATA  -0.033079, 
DATA  0.000063, 
DATA  0.002432, 


-0.007686, 

0.246045, 

-0.000417, 

-0.006519, 


-0.000053, 

0.007737, 

0.009682, 

-0.033848, 


-0.000209 

0.003644 

-0.004241 

0.126897 


FOR  L%  =  1  TO  4:  FOR  M%  =  1  TO  4 
READ  K#(L%,M%)  :  NEXT  M% 

NEXT  L% 


************************************************************ 
LOCATE  10,  20:  INPUT  "Type  the  name  of  the  voltage  output 
file";  VdL$ 

VOL$  =  "C:\LW\INSTR\DOUG\"  +  VOL$  +  ".DAT" 

OPEN  VOL$  FOR  APPEND  AS  #1 

CLS;  LOCATE  10,  20:  INPUT  "Type  the  name  of  the  forces  output 
file";  FM$ 

FM$  =  "C:\LW\INSTR\DOUG\"  +  FM$  +  ".DAT" 

OPEN  FM$  FOR  APPEND  AS  #2 

CLS:  LOCATE  10,  20:  INPUT  "TYPE  THE  NAME  OF  LIFT/DRAG  FILE"  ; 
LD$ 

LD$="C:  \LW\INSTR\DOUG\"+LD$+" .DAT" 

OPEN  LD$  FOR  APPEND  AS  #3 


500 

****************************************************** 
•ALPHA!  IS  AOA  READ  OFF  HUB  OF  TURNTABLE 
•ALPHA#  IS  AOA  OF  MODEL  WRT  AIRFLOW 

I  ***************************************************** 
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CLS:  LOCATE  10,  20:  INPUT  "Input  the  Test  AOA  (deg.)”;  alpha! 
ALPHA#*90  -  ALPHA! 

>  hit  it  1e-k1t1fk  It  Ifk -k  It  ielfk  It  It  Is -k  It  is-k  hit  *  hit -klele  It  Is -k  It  It -kit  ■kick -k  Is  *  It -k-k  is -kit  ***11 -k 

•CANARD  INCIDENCE  IS  WRT  TO  FUSELAGE  CENTERLINE 

*hhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhhh 


LOCATE  15,  20:  INPUT  "Input  canard  incidence  (deg.)";  CANARD# 

LOCATE  20,  20:  INPUT  "Input  oscillation  frequency  (Hz)";  HZ# 

CLS:  LOCATE  5,  20:  INPUT  "Is  this  a  tare  (zero  load)  reading? 
(Y/N)";  A$ 

IF  A$  »  "Y"  THEN  CALL  tare  (eanO#,eaaO#,ebnO#,ebaO#, alpha#) 

LOCATE  23,15:  INPUT  "Ready  to  take  readings?  (Y/N)";  B$ 

IF  B$  =  "Y"  THEN  CALL  volt  (ean#,eaa#,ebn#,eba#, alpha#) 

IF  B$  <>  "Y"  THEN  GOTO  5000 

•  Correcting  for  zero  load  values. 

eaa#  =  eaa#  -  eaaO# 
eba#  =  eba#  -  ebaO# 
ean#  =  ean#  -  eanO# 
ebn#  =  ebn#  -  ebnO# 

CALL  forces 

(K#  ( )  ,  eaa# ,  eba# ,  ean# ,  ebn# ,  AX# ,  Max# ,  NORM# ,  Mnorm# ,  alpha# ) 

CALL  aero  (AX#,NORM#, LIFT#, DRAG#, alpha#) 

PRINT  "  " 

PRINT  "  AOA  EAA  (mV)  EBA  (mV)  EAN  (mV)  EBN  (mV)" 

PRINT  ”  ***  hhhhhhhh  hhhhhhhh  hhhhhhhh  hhhhhhhh*^ 

PRINT  USING  "  ####.######";  alpha#;  eaa#;  eba#;  ean#;  ebn# 

PRINT  #1,  USING  "#####.######";  alpha#;  eaa#;  eba#;  ean#;  ebn# 

PRINT  "  " 

PRINT  "  AXIAL  (lb)  MOMax  (ft- lb)  NORMAL  (lb)  MOMnorm(ft-lb) " 
PRINT  "  hhhhhhhhhh  hhhhhhhhhhhhh  hhhhhhhhhhh  hhhhhhhhhhhh*^ 

PRINT  USING  "  ####.######";  AX#;  Max#;  NORM#;  Mnorm# 

PRINT  #2,  USING  "####.#####";  alpha#;  AX#;  NORM#;  LIFT#;  DRAG# 

PRINT  "  " 

PRINT  "  Lift  (lb)  Drag  (lb)  CANARD  (AOA)  CANARD  FREQ  (HZ)" 
PRINT  "  hhhhhhhhh  hhhhhhhhh  hhhhhhhhhhhh  hhhhhhhhhhhh** 
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PRINT  USING  '•  ####.######";  LIFT#;  DRAG#;  CANARD#;  HZ# 

PRINT  #3,  USING  "#####.######";  LIFT#;  DRAG#;  CANARD#;  HZ#; 
ALPHA# 

LOCATE  23,  15:  INPUT  "Do  you  want  another  reading?  (Y/N)”; 
ANS$ 

IF  ANS$  =  THEN  GOTO  500 

5000  CLOSE  #1 
CLOSE  #2 
CLOSE  #3 

END 

•  *********************************************************** 

•END  OF  MAIN  PROGRAM.  SUBROUTINES  START  FROM  HERE 
•  *********************************************************** 
************************************************************ 
•  *********************************************************** 

SUB  volt  (ean#,eaa#,ebn#,eba#, alpha#) 

•  *********************************************************** ' 

'  S/R  to  read  Channel  0,2, 4, 6  on  MIO-16L-9  for  Analog  Voltage 

t 


•  Setting  Board  code  for  MIO-16L-9 
board . code%»0 


************************************************************ 


errl.n\ua%  =  Init.DA.Brds(l,  board. code%) 


err2 .  n\im% 
err3 .  nu]n% 
err4.nuin% 
errs . num% 


AI. Setup (1,  0,  1) 
AI. Setup (1,  2,  1) 
AI.Setup(l,  4,  1) 
AI. Setup (1,  6,  1) 


'  Configure  and  set  clock  to  IMHZ 


err6.num%  =  CTR. Clock  (1,  1,  1,  1) 
err7.nuin%  =  CTR.Config  (1,  1,  0,  0,  0,  0) 


LWtotal!  «  0 


FOR  i%  »  1  TO  1000 


err8.num%  =  CTR.EvCount  (1,  1,  1,  0) 
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'  CH  0  =  Eaa 

err9.nuia%  =  AI.Read(l,  0,  1,  value0%) 

erl0.niua%  =  AI.Scale(l,  l,  valueo%,  eaa . array# ( i% ) ) 

•  CH  2  =  Eba 

erll.num%  =  AI.Read(l,  2,  1,  value2%) 

erl2.nuni%  =  AI. Scale (1,  l,  valued %,  eba. array# ( i%) ) 

•  CH  4  =  Ean 

erl3.nuin%  =  AI.Read(l,  4,  1,  value4%) 

erl4.num%  =  AI. Scale (1,  l,  value4%,  ean . array# ( i% ) ) 

'  CH  6  =  Ebn 

erl5.nim%  =  AI.Read(l,  6,  1,  value6%) 

erl6.niii&%  =  AI. Scale (1,  l,  value6%,  ebn. array# (i%) ) 

erl7.nuin%  =  CTR.EvRead  (1,  1,  overflo%,  tcount%) 

LWtotal!  =  LWtotal!  +  tcount% 

NEXT  i% 

CLS: LOCATE  5, 15: PRINT  "Total  Time  is  "  LWtotal J*lE-6" 
seconds . " 

CALL  Mean  (eaa.array#() ,  1000,  eaa#) 

CALL  Mean  ( eba . array# () ,  1000,  eba#) 

CALL  Mean  (ean.array#() ,  1000,  ean#) 

CALL  Mean  (ebn.array#() ,  1000,  ebn#) 


•  *********************************************************** 
'  This  multiplication  (*1000)  will  make  the  voltages  in  mV 

eaa#=eaa#*1000 

eba#*eba#*l000 

ean#=ean#*1000 

ebn#=ebn#*1000 

END  SUB 


I  *********************************************************** 
************************************************************ 


SUB  forces 

(K# ( ) , eaa# , eba# , ean# , ebn# , AX# , Max# , NORM# , Mnorm# , alpha#) 
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•  FORCES  AND  MOMENTS  CALCULATIONS  (See  thesis  for  explanation) 

AX#  =  K#(l,l)*eaa#  +K# ( 1 , 2 ) *eba#  +K# (1, 3) *ean#  +K#(1,4) *ebn# 
Max#  =  K#(2,l)*eaa#  +K# ( 2 , 2 ) *eba#  +K#(2,3) *ean#  +K#(2,4) *ebn# 
NORM#  =  K#(3,l)*eaa#  +K#(3,2)*eba#  +K# ( 3 , 3 )  *ean#  +K#(3,4)*ebn# 
Mnorm#  =  K#  ( 4 , 1 )  *eaa#+K#  (4,2)  *eba#  +K#  (4,3)  *ean#  +K#  (4,4)  *ebn# 

END  SUB 

************************************************************ 

************************************************************ 

SUB  aero  (AX#, NORM#, LIFT#, DRAG#, alpha#) 
************************************************************ 

PI#  =  3.14159265359 

LIFT#  =  AX#*COS(PI#/180*alpha#)  +  NORM#*SIN (PI#/180*alpha#) 
DRAG#  =  AX#*SIN(PI#/180*alpha#)  -  NORM#*COS (PI#/180*alpha#) 
END  SUB 

************************************************************ 
<  *********************************************************** 
SUB  tare  (ean#,eaa#,ebn#,eba#, alpha#) 

^  it  it  it  it  it  it  it  it  it  is  it  icit  it  it  it  icicieicitieic  it  It  it  ^  it  it  it  it  it  it  it  if  ie  it  it  it  it  icit  it  it  it  if  it  it  it  if  icit  it  it  ifkic 
t 

'  S/R  to  read  Channel  0,2, 4, 6  on  MI0-16L-9  for  Analog  Voltage 

I 

'  Setting  Board  code  for  MIO-16L-9 
board . code%=0 

•  *********************************************************** 

CLS:  LOCATE  5,  20:  INPUT  "Ready  to  take  tare  readings?  (Y/N)"; 
T$ 

IF  T$  <>  "Y"  THEN  RETURN 

errl.n\m%  =  Init.DA.Brds(l,  board. cod'' %) 

err2.n\UD%  «  AI. Setup (1,  0,  l) 
err3.nuin%  =  AI. Setup (l,  2,  1) 
err4.nuin%  *  AI. Setup (l,  4,  1) 
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err5.nuin%  =  AI. Setup (1,  6,  1) 

'  Configure  and  set  clock  to  IMHZ 

err6.nuin%  =  CTR. Clock  (1,  1,  1,  1) 
err7.nuin%  =  CTR.Config  (1,  1,  0,  0,  0,  0) 

LWtotal!  =  0 

FOR  i%  =  1  TO  1000 

err8.nvim%  =  CTR.EvCount  (1,  1,  1,  0) 

'  CH  0  =  Eaa 

err9.num%  =  AI.Read(l,  0,  1,  value0%) 

erl0.nuin%  =  AI. Scale (1,  1,  value0%,  eaa.array#(i%) ) 

'  CH  2  =  Eba 

erll.nvun%  =  AI.Read(l,  2,  1,  value2%) 

erl2.num%  =  AI. Scale (1,  1,  value2%,  eba. array# (i%) ) 

'  CH  4  =  Ean 

erl3.num%  =  AI.Read(l,  4,  1,  value4%) 

erl4.num%  =  AI.Scale(l,  1,  value4%,  ean. array# ( i%) ) 

'  CH  6  =  Ebn 

erl5.num%  =  AI.Read(l,  6,  1,  value6%) 

erl6.num%  =  AI.Scale(l,  1,  value6%,  ebn. array# ( i%) ) 

erl7.nuia%  =  CTR.EvRead  (1,  1,  overflo%,  tcount%) 

LWtotal!  =  LWtotal!  +  tcount% 

NEXT  i% 

CLSrLOCATE  5,15:PRINT  "Total  Time  is  "  LWtotal! *lE-6" 
seconds. " 

CALL  Mean  (eaa. array# () ,  1000,  eaa#) 

CALL  Mean  (eba.array#() ,  1000,  eba#} 

CALL  Mean  (ean.array#() ,  1000,  ean#) 

CALL  Mean  (ebn.array#() ,  1000,  ebn#) 

************************************************************ 

'  This  multiplication  (*1000)  will  make  the  voltages  in  mV 

eaa#=eaa# * 1 0 0 0 
eba#seba# * 1 0 0 0 
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ean#=ean#*1000 

ebn#=ebn#*1000 

PRINT  •'  '• 

PRINT  '•  AOA  EAA  (mV)  EBA  (mV)  EAN  (mV)  EBN(mV)" 

PRINT  "  ***  ********  ********  ********  *******" 

PRINT  USING  ••  ####.######";  alpha#;  eaa#;  eba#;  ean#;  ebn# 

PRINT  #1,  USING  "#####.######";  alpha#;  eaa#;  eba#;  ean#;  ebn# 

END  SUB 
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APPENDIX  D:  EXPERIMENTAL  DATA 
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MODEL  AOA= 


LIFT 


NEAM  LIFT  = 


STD.  DEV. 


3.182*STD.DEV 


LOWER  LIMIT 


FEAN  Cl  » 


Cl  STD.DEV.s 


3.182*STD.DEV.= 


UPPER  LIMIT*. 


LOWER  LIMIT* 


DRAG 


MEAN  DRAG  * 


STD.  DEV.* 


3.182*STD.DEV.= 


UPPER  LIMIT* 


LOWER  LIMIT* 


MEAN  Cd  * 


Cd  STD. DEV.* 


3.182*STD.DEV.* 


UPPER  LIMIT* 


LOWER  LIMIT* 


36.88612 

36.6991 

36.63613 

36.51826 

36.39341 

36.18825 

0.095674 

0.076712 

0.06343 

0.087532 

0.071285 

[« 

.074271 

0.304434 

0.244098 

0.201833 

0.278528 

0.22683 

0.236331 

37.19056 

36.94319 

36.83797 

36.79679 

36.62024 

36.42458 

36.58169 

36.455 

36.4343 

36.23973 

36.16658 

35.95192 

1.438407 

1.431113 

1.428658 

1.424061 

1.419193 

.411192 

0.003731 

0.002991 

0.002473 

0.003413 

0.00278 

0.002896 

0.011872 

0.009519 

0.007871 

0.010861 

0.008845 

0 

.009216 

1.450278 

1.440632 

1.436529 

1.434923 

1.428038 

1 

.420408 

1.426535 

1.421595 

1.420787 

1.4132 

1.410347 

1 

.401977 

14.62032 

14.56883 

14.5075 

14.46842 

14.30809 

0.036154 

0.041282 

0.027567 

0.03825 

0 

.021701 

0.032489 

0.115043 

0.13136 

.087718 

.121712 

0 

.069052 

0.10338 

14.92012 

14.75168 

14.65654 

14.62922 

14.53747 

14.41147 

14.69003 

14.48896 

14.48111 

14.38579 

14.39937 

14.20471 

0.577337 

0.570132 

B 

.568124 

B 

.565733 

.564209 

0.557957 

0.(K)141 

0.00161 

0.001075 

B 

.001492 

0.000846 

0 

.001267 

0.004486 

0.005122 

0.003421 

B 

.004746 

0.002693 

0 

.004031 

0.581823 

0.575255 

0.571545 

.570479 

0.566902 

0 

.561988 

0.572851 

0.56501 

0.564703 

0.560987 

0.561516 

0 

.553925 
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CANARD  AOA 


ANPLITUDE= 


DELTA  P  (cm.  H20)= 


q  <lbf/ft*2)= 


S  (ft*2)= 


MAC  (in.)- 


VELOCITY  <ft/s)= 


FREQUENCY 


Hz= 


LIFT 


MEAN  LIFT  » 


STD.  DEV.s 


3.182»STD.DEV.= 


UPPER  LIMITS 


LOUER  LIMITS 


MEAN  Cl  s 


Cl  STD.DEV.s 


3.182*STD.DEV.s 


UPPER  LIMITS. 


LOUER  LIMITS 


STD.  DEV.s 


3.182*STD.DEV. 


UPPER  LIMITS 


LOUER  LIMITS 


MEAN  Cd  s 


Cd  STD.DEV.s 


3.182*STD.DEV. 


UPPER  LIMITS 


LOUER  LIMITS 


5 


12 


27.33 


.9383 


5.38 


151.6454 


0.046439  0.092877  0.139316  0.185754  0.232193 


36.73228 


0.108599 


0.345563 


37.07784 


36.38672 


1.432407 


.004235 


0.013476 


1.445883 


1.418932 


37.22737 


37.24034 


37.12729 


36.95755 


.095732 

m 

.045358 

0 

.079247 

0.087175 

0 

0.304621 


37.53199 


36.92275 


1.451714 


.003733 


0.011879 


1.463593 


1.439835 


.144331 


37.38467 


37.09601 


1.45222 


0.001769 


0.005628 


1.457848 


1.446591 


0.252165 


37.37945 


36.87512 


1.447811 


0.00309 


0.009833 


1.457644 


1.437978 


.27739 


37.23494 


36.68016 


1.441192 


0.003399 


0.010817 


1.452009 


1.430375 


36.80912 


.06417 


.204187 


37.01331 


36.60493 


1.435404 


0.002502 


.007962 


1.443366 


1.427441 


15.29447 

13.97594 

13.95861 

13.93514 

13.8817  13.83226 

0.045459 

0.04185 

0.019704 

0.029195 

0.031175  0.01424 

0.14465 

0.133167 

0.0627 

0.092899 

0.099199  0.045312 

15.43912 

14.10911 

14.02131 

14.02804 

13.9809  13.87758 

15.14982 

13.84278 

13.89591 

13.84224 

13.78251  13.78695 

0.596421 

0.545004 

0.544328 

0.543413 

0.541329  0.539401 

0.001773 

0.001632 

0.000768 

0.001138 

0.001216  0.000555 

0.005641 

0.005193 

0.002445 

0.003623 

0.003868  U.001767 

0.602062 

0.550197 

0.546773 

0.547035 

0.545198  0.541168 

0.59078 

0.539811 

0.541883 

0.53979 

0.537461  0.537634 
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CANARD  AOA= 


AMPLITUDE 


DELTA  P  (on.  H20> 


q  (lbf/ft*2)= 


S  (ft*2>= 


MAC  (in.)= 


VELOCITY  <ft/s)= 


FREQUENCY 


Hz= 


LIFT 


•CAN  LIFT 


STD.  DEV. 


3.182*STD.0EV 


UPPER  LIMIT 


LOWER  LIMITS 


MEAN  Cl  s 


Cl  STO.DEV.s 


3.182*STD.DEV 


UPPER  LIMITS. 


LOWER  LIMITS 


MEAN  DRAG 


STD.  DEV. 


3.182*STD.DEV. 


UPPER  LIMIT 


LOWER  LIMITS 


MEAN  Cd  s 


Cd  STD.OEV.s 


3.182*ST0.DEV. 


UPPER  LIMITS 


LOWER  LIMITS 


0.046439  0.092877  0.139316  0.185754  0.232193 


36.62231 

36.36716 

36.29466 

36.22689 

36.07224 

35.88985 

0.098204 

0.10608 

0.062177 

0.02764 

0.079517 

0.070398 

0.312486 

0.337547 

0.197846 

0.087952 

0.253024 

0.224008 

36.93479 

36.70471 

36.49251 

36.31484 

36.32527 

36.11385 

36.30982 

36.02962 

36.09682 

36.13894 

35.81922 

35.66584 

1.428119 

1.418169 

1.415342 

1.412699 

1.406669 

1.399556 

0.00383 

0.004137 

0.002425 

0.001078 

0.003101 

0.002745 

0.012186 

0.013163 

0.007715 

0.00343 

0.009867 

0.008735 

1.440304 

1.431332 

1.423057 

1.416129 

1.416535 

1.408291 

1.415933 

1.405006 

1.407627 

1.409269 

1.396802 

1.390821 

15.49127 


0.034806 


0.110751 


15.60203 


15.38052 


0.604096 


.001357 


0.004319 


0.608415 


.599777 


15.2925 


.033264 


0.105847 


15.39835 


15.18666 


0.596344 


.001297 


.004128 


0.600472 


0.592217 


15.22074 


0.024964 


0.079437 


15.30017 


15.1413 


.593546 


.000974 


0.003098 


0.596644 


0.590448 


15.16547 


0.015401 


0.049006 


15.21447 


15.11646 


.591391 


0.000601 


0.001911 


0.593302 


0.58948 


15.15166 


0.035523 


0.113034 


15.2647 


15.03863 


.590852 


.001385 


.004408 


0.59526 


0.586444 


0.07283 


15.17518 


15.02952 


.588929 


0.000893 


.(K)284 


.5917 


0.586089 


80 


MODEL  ADA: 


LIFT 


MEAN  LIFT  ^ 


STD.  DEV.s 


3.182*STD.0EV.= 


LOUER  LIMITS 


MEAN  Cl  s 


Cl  STO.DEV. 


3.182*ST0.DEV.s 


UPPER  LIMITS. 


LOUER  LIMITS 


MEAN  DRAG 


STD.  DEV.s 


3.182*STD.0EV.s 


UPPER  LIMITS 


LOUER  LIMIT- 


MEAN  Cd  s 


Cd  STD.DEV.s 


3.182*STD.0EV.s 


UPPER  LIMITS 


LOUER  LIMITS 


5 

10 

15 

20 

.046439 

0.092877 

0.139316 

.185754 

0, 

25 


.232193 


36.88612 


0.095674 


35.62986 


0.114501 


35.46336 


0.081638 


35.28321 


0.076915 


.304434 

0 

.364342 

0.259771 

.244742 

0 

37.19056 


36.58169 


1.438407 


0.003731 


0.011872 


1.450278 


1.426535 


35.99421 


35.26552 


1.389418 


.004465 


0.014208 


1.403626 


1 .37521 


35.72314 


35.20359 


1.382925 


0.003184 


0.01013 


1.393055 


1.372795 


35.52795 


35.03847 


1.3759 


0.002999 


.009544 


1.385444 


1.366356 


35.09569 


0.121743 


.387385 


35.48308 


34.70831 


1.368587 


0.004747 


0.015106 


1.383694 


1.353481 


34.84926 


.133334 


0.424269 


35.27353 


34.42499 


1 .358977 


.005199 


0.016545 


1.375522 


1.342432 


14.80508 


0.036154 


0.115043 


14.92012 


14.69003 


6.577337 


0.(M)141 


0.004486 


0.581823 


0.572851 


14.65584 


.090177 


.286943 


14.94278 


14.3689 


0.571517 


0.003517 


0.01119 


.582707 


.560328 


14.49997 


14.42417 


14.34867 


.021084 

1  ° 

.01613 

0 

.05073 

II 

IE 

IE 


14.56705 


14.43288 


0.565439 


.000822 


0.002616 


.568055 


0.562823 


0.051324 


14.47549 


14.37284 


.562483 


.000629 


0.002001 


.564484 


.560482 


.161422 


14.51009 


14.18725 


.559539 


0.001978 


0.006295 


0.565834 


0.553244 


14.33285 


.047953 


0.152585 


14.48543 


14.18026 


.558922 


0.00187 


0.00595 


0.564872 


0.552972 
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MODEL  ADA 


LIFT 


MEAN  LIFT 


STD.  DEV. 


3.182*ST0.DEV 


UPPER  LIMIT= 


MEAN  Cl  = 


Cl  STO.DEV.- 


3.182*STD.0EV. 


UPPER  LIMIT*_ 


LOWER  LIMITS 


STD.  DEV.s 


3.182*STD.DEV. 


LOWER  LIMITS 


MEAN  Cd  s 


Cd  STD.DEV.s 


3.182*STD.DEV. 


UPPER  LIMITS 


LOWER  LIMITS 


36.73228 

35.88776 

36.04243 

35.80814 

35.8319 

35.53759 

0.108599 

0.049705 

0.080019 

0.18945 

0.090697 

0.115997 

0.345563 

0.158161 

0.254619 

0.602829 

0.288599 

0.369102 

37.07784 

36.04592 

36.29704 

36.41097 

36.1205 

35.90669 

36.38672 

35.7296 

35.78781 

35.20531 

35.5433 

35.16848 

1.432407 

1.399474 

1.405506 

1.39637 

1.397296 

1.385819 

0.004235 

0.001938 

0.00312 

0.007388 

0.003537 

0.004523 

0.013476 

0.006168 

0.009929 

0.023508 

0.011254 

0.014393 

1.445883 

1.405642 

1.415435 

1.419878 

1.40855 

1.400213 

1.418932 

1.393307 

1.395577 

1.372862 

1.386042 

1.371426 

15.29447 

15.19198 

15.07596 

14.9812 

14.95057 

14.84014 

0.045459 

0.029384 

0.02954 

0.063397 

0.034871 

0.050554 

0.14465 

0.093501 

0.093996 

0.201729 

.110958 

.160864 

15.43912 

15.28548 

15.16996 

15.18293 

15.06152 

15.001 

15.14982 

15.09848 

14.98196 

14.77947 

14.83961 

14.67927 

0.596421 

0.592425 

0.5879 

0.584205 

0.58301 

0.578704 

0.001773 

0.001146 

0.001152 

0.002472 

□ 

).00136 

■ 

.001971 

0.(K)5641 

0.003646 

0.003665 

0.007867 

0.004327 

0.006273 

0.602062 

0.596071 

0.591566 

0.592072 

0.587337 

0.584977 

0.59078 

0.588779 

0.584235 

0.576338 

.578683 

.572431 
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CANARD  AOAs 


M4Pl.ITUDE= 


DELTA  P  (cm.  H20) 


q  (lbf/ft^2)= 


S  <ft''2)= 


HAC  (in.}= 


VELOCITY  <ft/S>= 


FREQUENCY 


LIFT 


MEAN  LIFT  » 


STD.  OEV.» 


3.182*STD.DEV.- 


UPPER  LIMITS 


LOUER  LIMITS 


MEAN  Cl  s 


Cl  STD.DEV.s 


3.182*STD.DEV.s 


UPPER  LIMITS, 


LOUER  LIMITS 


DRAG 


MEAN  DRAG  s 


STD.  DEV. 


3.182*STD.DEV.s 


UPPER  LIMITS 


LOUER  LIMITS 


MEAN  Cd  s 


Cd  STD.DEV.s 


3.182*STD.DEV.s 


UPPER  LIMITS 


LOUER  LIMITS 


22 


10 


10 


12 


27.33 


.9383 


5.38 


151.6A54 


10 

15 

20 

.092877 

.139316 

.185754 

0 

25 


.232193 


36.62231 


0.098204 


34.80994 


0.139066 


35.27483 


0.069301 


35.19917 


0.06708 


35.04032 


0.022025 


.312486 

0.442508 

0.220516 

.213448 

IB 

.070084 

0 

36.93479 


36.30982 


1.428119 


35.25245 


34.36743 


1.357444 


35.49535 


35.05432 


1.375573 


.00383 

0.005423 

B 

.002702 

B 

0.012186 


1.U0304 


1 .415933 


0.017256 


1.3747 


1.340188 


0.008599 


1.384172 


1.366974 


35.41261 


34.98572 


1.372622 


.002616 


0.008324 


1.380946 


1.364299 


35.1104 


34.97023 


1.366428 


0.000859 


0.002733 


1.369161 


1.363695 


34.92356 


0.067171 


.213738 


35.1373 


34.70982 


1.361875 


0.002619 


0.008335 


1.37021 


35354 


15.39537 


0.157002 


15.10609 


0.063047 


15.08536 


0.025511 


.49958 

.200617 

0.081177 

B 

15.89495 


14.89579 


.006122 


0.019482 


.619837 


.580874 


15.30671 


14.90548 


0.589075 


0.002459 


15.16654 


15.00418 


0.588267 


0.000995 


15.09592 


0.024043 


.076506 


15.17242 


15.01941 


0.588679 


0.000938 


m 

IE 

IB 


.007823 

_ 1 

66 

0.002983 

.001791 

0 

0.596899 


.581252 


0.591432 


0.585101 


.591662 


.585695 


E 

B 

E 


15.05513 


0.014435 


0.045933 


15.10107 


15.0092 


0.587088 


.000563 


.588879 


0.585297 


15.01974 


.036104 


0.114883 


15.13462 


14.90485 


.585708 


0.001408 


.00448 


0.590188 


.581228 
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MODEL  AOA= 


AMPLITUDE= 


DELTA  P  (an.  H20): 


q  <lbf/ft''2)= 


S  (ft^2)= 


MAC  (in.)= 


VELOCITY  (ft/s)= 


FREQUENCY 


LIFT 


MEAN  LIFT  = 


STD.  DEV. 


3.182*STD.DEV 


UPPER  LIMIT 


LOWER  LIMIT- 


MEAN  Cl  - 


Cl  STO.OEV.s 


3.182*STD.DEV 


UPPER  LIMITS 


5 


12 


27.33 


.9383 


5.38 


151.6454 


5 

10 

15 

.046439 

0.092877 

.139316 

0 

.185754  0.232193 


41.7219 

41.96594 

42.09249 

41.72075 

41.57605 

41.58474 

0.255063 

0.093203 

0.140407 

0.105792 

0.102142 

0.146297 

0.811609 

0.296573 

0.446774 

0.33663 

0.325015 

0.465516 

42.53351 

42.26251 

42.53926 

42.05738 

41.90106 

42.05025 

40.9103 

41.66936 

41.64571 

41.38412 

41.25103 

41.11922 

1.626982 

1.636498 

1.641433 

1.626937 

1.621294 

1.621633 

0.009946 

0.003635 

0.005475 

0.004125 

0.003983 

0.005705 

0.031649 

0.011565 

0.017422 

0.013127 

0.012674 

0.018153 

1 .658631 

1.648064 

1.658856 

1.640064 

1.633969 

1.639786 

1 .595333 

1.624933 

1.624011 

1.61381 

mmmm 

1.60348 

MEAN  DRAG  s 

28.39872 

28.26335 

28.30194 

28.12068 

28.10789 

28.09581 

STD.  DEV.= 

0.108491 

0.033945 

B 

.066198 

0.05729 

0.052261 

0.052003 

3.182*STD.DEV.= 

0.345218 

0.108013 

B 

.210641 

0.182297 

0.166294 

0.165475 

UPPER  LIMITS 

28.74393 

28.37136 

28.51258 

28.30298 

28.27418 

28.26128 

LOWER  LIMITS 

28.0535 

28.15534 

28.0913 

27.93838 

27.9416 

27.93033 

MEAN  Cd  s 

1.107433 

1.102154 

1 

.103659 

1.09659 

1.096092 

1.095621 

Cd  STD.DEV.s 

0.004231 

0.001324 

0 

.002581 

0.002234 

0.002038 

0.002028 

3.182*STD.DEV.s 

0.013462 

0.004212 

B 

.008214 

0.007109 

0.006485 

0.006453 

UPPER  LIMITS 

1.120895 

1.106366 

1 

.111873 

1.103699 

1.102576 

1.102073 

LOWER  LIMITS 

1.093971 

1.097942 

1 

.095445 

1.089482 

1.089607 

1.089168 
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MODEL  AOA= 


CAMARO  AOA- 


LIFT 


MEAN  LIFT  = 


STD.  OEV.= 


3.182*STD,DEV 


UPPER  LIMITS 


LOWER  LIMITS 


MEAN  Cl  - 


Cl  STO.DEV.s 


3.182*STD.0EV 


UPPER  LIMITS 


LOWER  LIMITS 


DRAG 


MEAN  DRAG 


STD.  OEV.s 


3.182*STD.DEV 


UPPER  LIMITS 


LOWER  LIMITS 


MEAN  Cd  s 


Cd  STO.DEV.s 


3.182*ST0.DEV 


UPPER  LIMITS 


LOWER  LIMITS 


43.24928 


0.152075 


42.15297 


0.211123 


42.23697 


0.163877 


42.16644 


0.139068 


.483902 

0.671792 

55 

■ 

.442515 

43.73318 


42.76537 


1.686543 


0.00593 


.01887 


1.705413 


1.667673 


42.82476 


41 .48118 


1.643792 


0.008233 


0.026197 


1.669989 


1.617595 


42.75842 


41.71551 


1.647067 


0.006391 


0.020335 


1.667402 


1.626733 


42.60895 


41.72392 


1.644317 


.005423 


0.017256 


1.661573 


1.627061 


42.14693 


0.037684 


.119912 


42.26684 


42.02702 


1 .643556 


0.00147 


0.004676 


1.648232 


1.63888 


41 .85903 


.145572 


0.463211 


42.32224 


41.39582 


1.632329 


0.005677 


0.018063 


1 .650393 


1.614266 


28.14131 

27.77621 

27.7543 

0.026412 

0.105301 

0.092813 

0.084044 

0.335068 

0.295332 

28.22536 

28.11128 

28.04963 

28.05727 

27.44114 

27.45896 

1.097395 

1.083158 

1.082303 

0.00103 

0.004106 

0.003619 

0.003277 

0.013066 

0.011517 

1.100672 

1.096224 

1.09382 

1.094118 

1.070091 

1.070786 

27.69721 


.276451 


27.97366 


27.42076 


1.080077 


.003388 


.01078 


1.090857 


1.069296 


27.71591 


.007353 


.023398 


27.73931 


27.69251 


0.000287 


0.000912 


1.081718 


1.079894 


27.5449 


.072445 


.230519 


27.77542 


27.31438 


1.074137 


0.002825 


0.008989 


1.083127 


1.065148 
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MODEL  AOA= 


CANARD  AOA= 


DELTA  P  (on.  H20) 


q  (lbf/ft'2)= 


S  (ft''2)= 


MAC  (in.)- 


VELOCITY  (ft/S)= 


FREQUENCY 


MEAN  DRAG  > 


STD.  DEV.s 


3.182*STD.DEV.‘ 


UPPER  LIMIT 


LOUER  LIMIT 


MEAN  Cd  « 


Cd  STD.DEV.s 


3.182*STD.DEV.> 


UPPER  LIMIT* 


LOUER  LIMIT* 


34 


-10 


5 


12 


27.33 


.9383 


5.38 


151 .6454 


10 

15 

20 

.092877 

.139316 

0.185754 

0. 

25 


.232193 


MEAN  LIFT  * 

42.40407 

41.79588 

41.92224 

41.78692 

41.88277 

41.62126 

STD.  DEV.* 

0.059618 

0.136226 

0.160787 

0.158283 

0.055838 

0.023508 

3.182*STD.DEV.= 

0.189706 

0.43347 

0.511623 

0.503657 

0.177678 

0.074802 

UPPER  LIMIT* 

42.59378 

42.22935 

42.43386 

42.29057 

42.06044 

41.69606 

LOUER  LIMIT* 

42.21437 

41.36241 

41.41062 

41.28326 

41.70509 

41.54646 

MEAN  Cl  * 

1 .653584 

1.629867 

1.634794 

1.629517 

1.633255 

1 .623057 

Cl  STD.DEV.* 

0.002325 

0.005312 

0.00627 

0.006172 

0.002177 

0.000917 

3.182*ST0.0EV.* 

0.007398 

0.016904 

0.019951 

0.019641 

0.006929 

0.002917 

UPPER  LIMIT*, 

1.660982 

1.64677 

1.654746 

1.649158 

1.640184 

1 .625974 

LOUER  LIMIT* 

1.646186 

1.612963 

1.614843 

1.609877 

1.626326 

1 .620141 

27.32982 

26.89041 

26.93577 

26.84227 

26.93257 

26.75407 

0.033418 

0.075427 

0.074768 

0.091834 

0.02875 

r? 

.015526 

0.106336 

.240008 

B 

.237913 

0.292216 

0.091482 

0.049403 

27.43616 

27.13041 

27.17368 

27.13449 

27.02405 

26.80348 

27.22348 

26.6504 

26.69786 

26.55005 

26.84109 

26.70467 

1.06575 

Q 

.048615 

.050384 

1.046738 

1.050259 

Q 

.043298 

0.001303 

B 

.002941 

B 

.002916 

0.003581 

0.001121 

0.000605 

0.004147 

0.009359 

0 

.009278 

0.011395 

0.003567 

0 

.001927 

1.069897 

1 

.057974 

1 

.059661 

1 .058133 

1.053827 

1 

.045225 

1.061604 

1 

.039256 

1 

.041106 

1 .035343 

1.046692 

1 

.041372 
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MODEL  AOA= 


AMPIJTUDE= 


DELTA  P  (cm.  H20) 


q  (lbf/ft"2)= 


S  (ft*2)= 


MAC  (in.)= 


VELOCITY  (ft/s)= 


FREQUENCY 


LIFT 


MEAN  LIFT  = 


STD.  DEV. 


3.182*STD.DEV. 


UPPER  LIMITS 


MEAN  Cl  3 


Cl  STO.DEV.s 


3.182*STD.DEV. 


UPPER  LIMIT= 


LOWER  LIMITS 


DRAG 


MEAN  DRAG 


STD.  DEV. 


3.182*STD.DEV.= 


UPPER  LIMITS 


LOWER  LIMIT 


MEAN  Cd  s 


Cd  STO.DEV.s 


3.182*STD.DEV.s 


UPPER  LIMITS 


LOWER  LIMITS 


41.7219 

41 .68476 

41.61744 

41.4616 

41.51339 

41 .33265 

0.255063 

0.045258 

0.119078 

0.037197 

0.117325 

0 

.127482 

0.811609 

0.144012 

0.378907 

0.11836 

0.373328 

0 

.405648 

42.53351 

41.82877 

41.99635 

41.57996 

41.88671 

41.7383 

40.9103 

41.54074 

41.23854 

41.34324 

41.14006 

40.927 

1.626982 

1.625533 

1.622909 

1.616831 

1.618851 

1 

.611803 

0.009946 

0.001765 

0.004644 

0.001451 

0.004575 

0 

.004971 

0.031649 

0.005616 

0.014776 

0.004616 

0.014558 

0.015819 

1.658631 

1.631149 

1.637684 

1.621447 

1.633409 

1 

.627621 

1.595333 

1.619918 

1.608133 

1.612216 

1.604292 

1 

.595984 

28.39872 


.108491 


0.345218 


28.74393 


28.0535 


1.107433 


.004231 


.013462 


1.120895 


1.093971 


28.09699 


0.026359 


0.083873 


28.18087 


28.01312 


1.095667 


0.001028 


0.003271 


1 .098937 


1.092396 


27.85282 


.051637 


.16431 


28.01713 


27.68851 


1.086145 


.002014 


0.006407 


1.092553 


1.079738 


27.76555 


.027289 


.086833 


27.85239 


27.67872 


1.082742 


.001064 


0.003386 


1.086128 


1.079356 


27.74751 


.053095 


0.168949 


27.91646 


27.57856 


1.082038 


0.00207 


0.006588 


1.088627 


1.07545 


27.67066 


0.062103 


.197612 


27.86827 


27.47305 


1.079042 


0.002422 


0.007706 


1.0867 


1.071336 


MODEL  AOA 


DELTA  P  (cm.  H20}= 


q  (lbf/ft"2)= 


S  (ft*2)= 


MAC  (in.)= 


VELOCITY  (ft/s)= 


FREQUENCY 


10 


12 


27.33 


0.9383 


5.38 


151.6454 


0.046439  0.092877 


15 


.139316 


.185754  0.232193 


LIFT 


MEAN  LIFT 


STD.  DEV. 


3.182*STD.DEV. 


UPPER  LIMIT 


LOUER  LIMIT 


Cl  STO.OEV.s 


3.182*STD.DEV. 


UPPER  LIMITS 


LOUER  LIMITS 


43.25268 

41.06015 

41.09099 

41.15339 

41.01945 

41.12252 

0.14652 

0.10197 

0.093618 

0.151984 

0.072549 

0.161978 

0.466227 

0.324469 

0.297892 

0.483613 

0.230852 

0.515414 

43.71891 

41.38462 

41.38888 

41.637 

41.25031 

41 .63794 

42.78645 

40.73568 

40.7931 

40.66977 

40.7886 

40.60711 

1.686676 

1.601176 

1.602379 

1.604812 

1.599589 

0.005714 

0.003976 

0.003651 

0.005927 

0.002829 

0.006316 

0.018181 

0.012653 

0.011617 

0.018859 

0.009002 

0.020099 

1.704857 

1.613829 

1.613996 

1.623671 

1.608592 

1.623708 

1.668495 

1.588523 

1.590762 

1.585953 

1.590587 

1.58351 

DRAG 


MEAN  DRAG  = 


3.182*STD.DEV. 


UPPER  LIMITS 


LOUER  LIMITS 


MEAN  Cd  s 


Cd  STO.DEV.s 


3.182*STD.DEV. 


UPPER  LIMITS 


LOUER  LIMITS 


28.14131 


.026412 


.084044 


28.22536 


28.05727 


0.00103 


0.003277 


1.100672 


1.094118 


27.4693 


0.059916 


.190654 


27.65995 


27.27864 


1.071189 


0.002336 


.007435 


1.078624 


1.063755 


27.24352 


27.30029 


27.29642 


.064487 

.077418 

0 

.044499 

0.205199 


27.44872 


27.03833 


1.062385 


0.002515 


.008002 


1.070387 


1.054383 


0.246346 


27.54664 


27.05395 


1.064599 


.003019 


.009606 


1 .074205 


1 .054992 


0.141597 


27.43802 


27.15482 


1.064448 


.001735 


0.005522 


1.069969 


1.058926 


27.31911 


.083987 


0.267247 


27.58636 


27.05186 


1.065333 


0.003275 


0.010422 


1.075754 


1.054911 
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MODEL  AOA= 


34 


CANARD  AOA= 


AMP-.ITUDE 


DELTA  P  (cm.  H20) 


q  (lbf/ft"2)= 


S  (ft^2)= 


MAC  (in.)s 


VELOCITY  (ft/s)= 


FREQUENCY 


Hz= 


LIFT 


MEAN  LIFT  = 


3.182*STD.DEV. 


UPPER  LIMIT= 


LOWER  LIMIT= 


MEAN  Cl  = 


Cl  STD. DEV. = 


3.182*STD.DEV. 


UPPER  LIMIT= 


LOWER  LIMITS 


DRAG 


MEAN  DRAG  - 


STD.  DEV.= 


3.182*STD.DEV 


UPPER  LillT 


LOWER  LIMITS 


MEAN  Cd  s 


Cd  STD. DEV. s 


3.182*STD.DEV 


UPPER  LIMITS 


42.40407 


0.059618 


0.189706 


42.59j,o 


42.21437 


1.653584 


0.002325 


0.007398 


41.63824 


41.53664 


41.77239 


.15827 

■ 

.120164 

.061498 

r° 

1.646186 


.503615 


42.14185 


41.13462 


1.623719 


0.006172 


0.019639 


1.643358 


1. 


0.382363 


41.919 


41.15428 


1.619758 


0.004686 


0.014911 


1.634668 


1.604847 


0.195687 


41. 


41.57671 


1.628951 


0.002398 


0.007631 


1.636582 


1.62132 


41.56922 


.082291 


0.261851 


41.83107 


41 .30736 


1.621028 


0.003209 


0.010211 


1.631239 


1.610817 


27.32982 


0.033418 


.106336 


27.43616 


27.22348 


1.06575 


.001303 


.004147 


1.069897 


1.061604 


27.57928 


0.10223 


0.325295 


27.90457 


27.25398 


1.075478 


0.003987 


0.012685 


1.088163 


1.062793 


27.31677 


0.064667 


0.20577 


27.52254 


27.111 


1.065241 


.002522 


0.008024 


1.073266 


1.057217 


27.45698 


.031294 


.099577 


27.55655 


27.3574 


1.070709 


0.00122 


0.003883 


1.074592 


1.066826 


27.4311 


0.048506 


0.154346 


27.58545 


27.27675 


1.0697 


0.001892 


.006019 


1.075719 


1.063681 


41 .6576 


0.093878 


0.298721 


41.95632 


41 .35888 


1 .624474 


0.003661 


0.011649 


1.636123 


1 .612826 


27.45643 


0.043012 


0.136864 


27.5933 


27.31957 


1.070688 


0.001677 


0.005337 


1.076025 


1 .06535 
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